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PART  II 

ASYMPTOTIG  EXPANSIONS  OF  SOLUTIONS 
OF  DIFFERINTIAL  EQUATIONS 


GHAPTER  1 

SOLUTION  OF  tm  RAOIAL  DiFFroENTIAL  EQUATION  FOR  TH  WAVES 

GONVERGENT  SERIES  SOLUTIONS  AROUND  x  s  0 
We  Gbnslder  equation  (1*^50),  PART  It 

*  ISS'A) I'Sgf  y  *  xfs  *  ^^3  R<*)  =  0  •(!  l-SO) 

The  origin  x^O  is  a  regular  singularity  of  this  equation.  In  order 
to  obtain  convergent  series  solutions  around  this  point,  we  follow 
the  classical  method  of  Frobenius  (9  pp*  396o*404).  First,  write 
the  equation  in  the  formt 

ix-<f (a+b)x+ab]x-  r*(x)+cx-  Rl[x)'f[x^+2ax^+(a^*v(r+l) )x^-v(v+l) • 
t (a*b)x=v(v+l)ab]  R(x)  =  0  • 

Substituting  the  formal  series  R(x)  =  equation 

yields:  a^jX^^^f (x,n+o)  =  0,  where: 

1-1 


1-2 


f  ( X , ii+cr) = ( n+t)  { a^fe) x+ala] + ( n+er)  cx+x^+2ax^+C a^-v ( v+1 )  ]x^- 

-v(v*fl)  (a+b)x-v(v+l)ab  =  fg(n+tr)+fj^(ii'*‘or)x+fg(n+b)x^+fj(n+df)x^+ 
+f^(n+Gf)x^* 

The  fimctlons  appeas'lng  In  last  expresslbn  are  defined  as 


follows i 

f^(n+o)  =  abC(n+a)  (n+ai‘l)-v(v+l)]=ab(n+ff*fv)  (n+a-v-1)  (1-1) 

fj^(n+a)  =  (a+b)  (n+e)  (n+o'-l)+e(n+e)-(a+b)v(v4‘l)  (1-2) 

fg(n+o)  =  (n+e)(n+e-l)+a^-v(yel)  (1-3) 

fj(a+e)  =  2a  (1-4) 

f^(n+cr)  =  1  .  (1^5) 

Indlclal  equatioa* 

f^(q')  =  ab(qf+v)(e-v-»l)  =  0  (1-6) 

with  roots: 

V+l  ,  CTgS  -V  .  (1-7) 

Reccurence  formula  for  the  coefficients; 


a^f  p  ima)  ^  (n+o-l )  +aj^^2f  g  (n+a-2)  +an»3^5  (n+cr-3 )  + 

^Sjj^^f ^(n+0«4)  s  0  5  0»ffi^l»2t3f»*»*  •  (1—8) 

V  is  a  positive  parameter»  since  negative  values  of  v  do  not  yield 
modes  independent  from  the  ones  corresponding  to  positive  values 
of  V*  We  have: 

§2  -  2v#l 


(1-9) 
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If  av+l  is  not  an  Integei',  the  method  yields  two  independent 
solutions.  In  pantieular*  we  define  as  H^(x)  and  Rg(x)  the 
solutions  Gorrespondlng  to  0^^=  v+l  and  -v*  respectively, 
with  ag=  1  ,  i.e. 

(x)  =  a^x^]  ,  :|x|  ^  inin(  |al  ,  |b|  )  (I-IO) 

pO 

Rg(x)  =  X  '^[1+^^  *  [xj  ^  min(  l-aj  »  (fe|  )  ,  2y*l^tnteger  (li44) — 

POX’  the  coefficients  a^^  we  use  (1-8)  with  orsc^sv+l  and  a  =  1, 
fox*  the  bj^'s,  again  (1»8)  with  osCgS-v  and  1*  We  can  find 
that: 

Alao: 

min(|a],|bl)  -  t  |^j  jt)<ja|  •  (1-12) 

We  observe  that  R|(0)  =  0  and  R2(0)  =  oo  ,  x=o  is  a  branch  point 
for*  both  solutions. 

If  2v^l  is  equal  to  a  positive  integer,  (1-10),  together 
with  (1*8),  continues  to  define  R^(x).  The  second  independent 
spLutlon  Rg(x)  in  this  case,  has  a  logarithmic  singularity  at 
X5p  and  is  going  to  be  found  later. 

Analytic  Continuation  of  R^(x) .  R-(x)  in  the  Right-Half 
x-Plane;  We  obtained  two  series  representations  for  the  solutions 
Rj^(x)  and  R2(x)  of  equation  (I  1-50).  They  are  valid  only  within 
the  circle  ]xl<mln(|al  ,  |bl  ),  where  they  both  converge  uniformly 
Md  absolutely.  We  proceed  to  obtain  representations  for  these 
PM*tleular  fvmctions  valid  in  the  whole  right-half  plane,  where 
X  ’varies  from  0  to  po  . 

We  use  a  bilinear  transformation  of  the  independent  variable 
in  the  form: 


1-4 


t 


x+p 


(1-13) 


The  eohStaht  p  is  cohvenlehtly  chosen  in  each  ease  so  as  to 
optimize  the  rate  of  convergence  of  the  resulting  series  in  t» 

The  conformal  mapping  of  x-plane  on  to  t-plane>  in  accordance 
with  (1-13) »  is  showi  in  figure  (1-1).  The  half-plane  to  the 
right  of  the  perpendicular  bisecting  the  line  between  o  and  -p 
(shaded  in  figure  (l«l))y  maps  within  the  tmit  circle  Itl^  1  of 
the  t-plane.  in  this  plane  x  varies  from  0  to  esa  along  a  straigdit 
line  in  the  fourth  quadrant,  or  along  the  real  axis.  Corresponding 
ly,  t  varies  from  0  to  1  within  the  imit  circle.  The  singular 
points  of  equation  (I  1-50)  map  as  follows t 
X  =  0  t  =  0 

X  s  -a  t  =  a/(a-p) 

X  =  -b  t  m  b/(b-p) 

X  =  o«  t  a  1  • 

The  equation  becomes: 


d^R  r  P®  ^ 

dt^  L  t ( P-a ) t+a] i ( p-b ) t+b i  t-lj 


dt  * 


r 

(p-a)t+a  p‘ 


1— 


h)t+b  (t*l)^  t^(t-l)^J 


R  =  0 


(1-14) 


It  now  has  three  regular  singularities  at  t=G,  t=b/(b-p), 
t=a/(a»p)  and  an  Irregular  singularity  at  t=l.  Thus,  the  trans- 
forma  tlon  has  preserved  the  nature  and  rank  of  the  singularities > 
mapping  them  at  their  image  points  in  the  t-plane*  This  Is  a 
well-teiowi  property  of  bilinear  transformations  (9  p.  437).  If  p 
is  chosen  as  shown  in  figure  (l-l) ,  the  singularities  x=9a  i^d 
x=-b  map  outside  the  unit  circle  in  the  t-plane.  A  power  series 
solution  of  equation  (1-14)  arotoxd  t=P  will  converge  uniformly 
and  absolutoly  for  |t|<l  and  would  provide  a  representation  valid 
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over  the  whole  laterval  of  Interest  O^Hex^co.  The  nuaerical 
computations  of  chapter  3*  PART  l,  showed  that  such  series  can  be 
used  for  an  ixL  3  or  4  times  larger  than  the  |x;l'  of  (l-lo) 
^d  (1-11) i  It  is  not  necessary  to  map  x=«a  and  x@«b  outside  the 
Circle  lt|  =  1»  Depending  on  the  case  under  consideration,  a  value 
of  p  can  be  chosen,  that  maps  either  or  both  of  these  points 
inside  |t|  =  hut  which  yields  better  series  in  t.  The  radius 

of  convergence  is  news  jtl  =  ||fp)  »| flp)  1  5  nevertheless, 

larger  values  of  x  can  be  used  with  such  series,  that  may  carry 
one  farther.  In  this  connection,  notice  that  x  varies  only  along 
a  straight  line  in  the  fourth  quadrant,  or  along  the  real  axis, 
from  0  to  60  and  that  one  is  only  interested  in  reaching  up  to 
such  values  of  x  after  which  the  asymptotic  expansions  of  Rj^(x) 
and  can  be  used.  The  fact  that  a  series  expansion  can  be 

found,  which  is  valid  for  all  x  from  0  to  e>o  ,  has  no  practical 
Significance  from  the  computational  point  of  view*  A  series  valid 
within  a  smaller,  finite  region  may  be  better  if  it  can  be  used 
for  larger  values  of  x.  These  facts  were  verified  numerically. 

|P|  should  be  chosen  as  large  as  possible,  but  not  as  large  as  to 
result  in  a  very  small  convergence  circle  in  the  t^plane.  In  each 
particular  case  an  optimvm  value  of  p  exists. 

For  all  Cases  t  to  VI,  Chapter  3,  PART  I,  p  was  given  the 
following  value: 

.  *  =  f!|  •  Cl-15) 

We  will  develop  the  series  expansions  In  t  for  %(^)  and  RgCx) 
for  this  particular  choice.  For  a  different  P  the  analysis 
follows  identical  lines.  Now,  x=..a  maps  on  t=^l  and  (1«14) 
becomes • 


•t%H)+t|||f(t+l)^tW(v+l)  (t+i)(t+||^)]R(t>  =0  .  (1-16) 

po 

Inset^tlftg:  R(t)  =  ^  equation  becomeit 

n=o  ” 

CMO 

%  (t,n+o)  =  0  ,  where  s 

n=6 

f  ( t  j  n+cr)  =  (a+o)  (a+o-1 )  C  t^+  ( h-5 )  t® + ( 9-3h)  t'^+g  ( h+i )  ( 2h-3 ) 

+  (l-3h)  tH^hl-f  (n+of){  2t^+  (h-3)  t^“4t^+(lG-6h)  t^+  ( 8h-6)  t^-(3h-l )  t]+ 

+{  ga^  ( h+1 )  -V  ( v+1 )  ]  t^+{  4a^  ( h+1 )  ( 1-h)  v  ( v+1 )  ]  t^+  ( 1+h)  £  20^+ v  ( v+1 )  ] 
+  (h-l)v(v+l)  t»hv(v+l)  =  2  fg^(a+cr)t®  . 

The  following  deflnltlona  were  uaedt 


h  =  b/(2a*b) 

(1-17) 

f^(x)  =  h[x(3 

c-l)-v(v+l)i  m  h(x+v)(x-v-l) 

(1-18) 

t^(x)  =  (lOl 

l)x^+(h-l)v(v+l) 

(1-19) 

fgCx)  =  (2h-^ 

5 ) x^+ ( 6h-3 ) x+ ( h+1 ) ( 2a^+ v^+ V ) 

(1-20) 

f^(x)  =  2(h+] 

L)x^+a(l-h)x+4a^(h+l)+(l-h)v(v+l) 

(1-21) 

tj^(x)  =  (2-31 

i)  (x^-x)-4x+2a^(h+l)-v(v+l) 

(1-22) 

fgCx)  =  (h-3) 

x^ 

(1-23 ) 

P 

fg(x)  S  X-+X 

(1-24) 

Indie lal  equ{ 

itlon*  ^0  with  roots: 

=  v+i 

,  cTg  =  -y 

(1-25) 

Reeurrenee  formulae 


i 
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Z  ^  ®  *  ®6~  ^  ‘  ®  ^  ^  1»2,^..  (1-26) 


m=o 


For  9  =  9|  =  v-t^l  we  always  obtalzi  a  solution t 


R(t)  =  Kt®[l+  2  s^t”3  , 

n=l  ^ 


(1-27) 


where  K  is  a  Gonstaat  and  the  eoefflclents  e„  are  determined 

n 

using  (1-26)  with  9  =  v+1.  If  9^-9g=  2v+l  is  hot  an  integer,  use 
of  (1-26)  and  (1-27)  with  9  =  -v  yields  a  second  ladependent 
solution.  If  2v‘fl  is  equal  to  a  positive  integer,  the  second 
solution  is  logarithmic  and  will  be  found  later. 

For  n=l  equation  (1-26)  yields:  e^=  -f2^(9)/fg(9+l) ;  since  in 
any  case:  9(9-1)  =  v(v-fl),  we  obtain: 


®1®  9^ 


1„ 

3h 


(1-28) 


As  x—0  and  t-»Oj  we  have: 


^  ^  X42a  as  |+(x72a)  "  ♦••••J 

Substituting  in  (1-27)  we  obtain: 


R(t)®K(|-)^(i.|.)*ni*(9*  i-Ij)  fj  (i.0(^,)*0(^2,] 


-  K 


i  •  k% 


,  K 


K 


x^tl-  ||p!+0(x8,] 


facing  K  =  (2a)®  and  referring  to  (1-10) , (1-11)  and  the  remarks 
following  them,  we  see  that  the  two  solutions  defined  by  (1-27), 
(1-26)  can  be  identified  with  Rj^(x)  and  ^^(x)  of  (1-10),  (1-11) 
In  the  corresponding  eases  9^-  V4l  and  -v  .  That  is: 


1-8 
(1-29) 

Solution  R|  Is  given  fon  o^v+l,  for  Os-v,  if  iv+l  is  not  an 
integer.  If  2v+l  is  equal  to  an  integer j  Rg  will  ne  found  later * 
The  series  in  x  eonverge  uniformly  and  ahsolutely  for 
|x|<min(  lal  ,  Ihl  ),  those  in  t  for  |t|<min(l,  lh|  )*  providing  the 
walytie  eontinuation  of  R^(x),  R^Cx)  into  the  right  half -plane 
Shown  in  figin*e  (i-l). 


Rj^de)  =  = 

g  n=l 


of^qri 


t^Cl+1  e^t"] 
n=l  ^ 


2v4l  IS  EQUAL  TO  AN  INTEOER 

In  this  case  the  second  independent  solution  R2(x)  of 
equation  (I  1-50)  becomes  logarithmic.  The  preceding  analysis 
is  validj  without  alteration*  to  the  solution  R2(x)  with  o^v+l* 

In  the  biconical  antenna  of  PART  I,  as  well  as  in  other  problems* 
the  second  solution  R2(x)  appears  indirectly.  What  is  actually 
directly  involved  is  R^(x) *  the  solution  which  satisfies  the 
radiation  condition  at  infinity*  corresponding  to  an  Integral 
value  of  V.  Evaluation  of  R^(x)  for  small  x*  where  its  asymptotic 
series  representation  can  not  be  used,  involves  both  R^(x)  and 
Rg(x)  of  integral  order.  An  equation  like  (1-54),  PART  I,  must  be 
available  to  provide  the  analytic  continuation  of  R^(x)  in  the 
vicinity  of  x=0.  Thus,  not  only  RgCx),  but  its  complete  asymptotic 
expansion  for  large  |x| ,  equation  (1-53),  PART  I,  must  be  known. 

Following  Whittaker  and  Watson  (16  pp.  200-201),  we  put: 

R2(x)  «  inx  Ri(x)  +  S(x)  ,  (1-50) 

oo 

where  R^^lx)  =  2  1,  is  the  first  solution  obtained 

~  m=o 

in  the  preceding  section*  and: 


U9 

S(x)  =  % 
n=d 

We  write  equation  (I  1«50)  in  the  following  formt 
f^(x)x%'‘(x)ft^(x)xR^(x)+‘f2(x)R(x)  =  0  *  (1^32) 

where : 


iLq(x)  =  x^»(a4b)x»ab 
t*.  (x)  =  cx 


f^(x)  =  x^+2ax^+[a^»v(v+l) 3x^« 
Substituting  (1*30)  into  (1*32) 


f^(x)x^CH^(x)lnx+2  ^ 


li.'(x)  R.  (x) 


v  (v‘*>l)  (a‘»-b)x*v(v+l)ab  * 

we  obtain: 

•  ^S^'(x)  ]H‘r^(x)x[R^(x)lnx+ 


R,  (x) 

♦  +s1ix)]+t2(^)C^i(^)3.nx+s(x)i  =  0  , 

or*  since  R^(x)  satisfies  (1*32): 


far*33) 

(1*34) 

(1*35) 


T^(x)  X%  "(  x)  +t^  (  x)  xS'(  x)  ♦Tg  (x)  S  (  x)  =5 

-  T^(x)[R^(x)-2xR^(x)]*‘f2^(x)R^(x)  i  (1*36) 

pp 

Substituting  (1*31)  ,(1*33) -(1*35)  and  Rj^(x)=  Z  in  the 

above  equation,  we  obtain: 


Z  b_x®"''f  (x,n*v)  =  £  a_x®-"''*^(x.in+v+l)  ,  (1*37) 

n=o  “  m=o  ® 

where : 

f(x,n*v)  =  (n*v)  (n*v*l)Cx-+(a+b)x'fab3+(n*v)cx+x\2ax^+[a^* 

O  ^  m 

-v(v-fl)]x^*v(v+l)(a+b)x*v(r4l)ab  =  Z 

m=P  - 


and: 
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P(x,ffl+v%l)  =  -C2(m+v+l)-l3Ex^%(a+b)x%ab]*cx  =  ^  P.Cm+v+Dx-^  . 


J=6  J 

fna  following  deflnitlonB  were  used: 

fg(n«v)  =  abC (n-v) (n-v-l)»v(v+i)l  =  abn(n«2v-l)  (1-38) 

f|(n-v)  =  (a+b) (n«v) (n-v-i)+(n-v)c-(a+b)v(v*l)  (1*39) 

fg(n*v)  =  ( n- V )  ( n*  V  *1 )  +a^-v  ( v^+il  (1*40) 

■l 

t^(n^v)  s  2a  (1*41) 

f^(n*v)  =  1  (1*42) 

P^(a+v+l)  =  -ab[2(m+v+l)*l]  (1*43) 

P^(m+v+l)  s  *(a+b)[2(m+v+l)*l]*c  =  *2(a+b) (m+v4l)+2b  (1*44) 

Pg(m+v+l)  =  *2(m+v+l)+l  .  (1*45) 

iquation  (1*37)  can  be  written: 

^  fen^f(x»n*v)  s  ajjjX®p(x»a+v+l) 


nso  -  aeo 

and  shows,  Inc Inden tally,  that,  with  f(x,n*v)  and  F(x,m+v4»l) 
being  pol^oalals  In  x.  It  can  not  be  satisfied  xuiless  2v^l  is 
equal  to  a  positive  Integer.  In  its  right-hand  side  put: 

B+2V+1  m  n  •  (1*46) 

It  becoaes: 

•o  ^ 

E  b  x-f(x,n«v)  =  S  ix^  p(x,n«v)  , 

n=o  “  n=2v+i 

Since  a^^=  a_2=.....=  0  ,  we  can  also  write  It  as  follows: 


1*11 

(x,n«v)  =  an-2v*l^*^  *  (1*^7) 

ns©  n— 6 

The  lowest  exponent  of  x  is  n=0  and  Its  eoefflcient  is  b^fg(0*v)= 
=0.  Generally  for  l^n^gv  we  have: 

Vo  *\-l^l  ^  ^  +\*2^2  ^  ^  ^ 

♦h^_^f^(n«v-4)  =  0  ;  0  ,  m  =  1,^,3m««  *  (1*48) 

For  all  these  values  of  n  we  observe  that  f^(n*v)  ^  0*  permitting 

to  use  successively  the  recurrence  formula  and  obtain  b^ «  b*> 

Pv+1'*'  ^ 

4i.i  bgyi  For  n=2v>^i,  i.6i  for  the  coefficient  of  x  in 
(1*47) »  we  have: 

=  aoFQ(v+l)  =  *ab(2vH^l)  .  (1*49) 

NOW,  however,  fQ(v+l)  ^  0,  according  to  (1*38).  In  order  to 
satisfy  (1*49),  vfe  choose  conveniently  the  value  of  bp*  left 
vindetermined  so  far*  Notice  that  b^,  bg^,  determined  so 

far  with  the  use  of  (1*4$) »  are  all  proportional  to  b^.  in  fact, 
we  can  write: 


Wn  •  V‘  ■  . 

and  use  (1*48),  with  d^^  in  place  of  b^,  to  determine  the  numbers 
do»  d^,  dg,....,  dgy  completely,  with  initial  conditions:  4^=  1, 
d^j^s?  d_g=.. ..=  0.  Then,  (1*49)  is  satisfied  if  we  take: 


ab(2v+l) 

a2,fj(v)+dj^.jfj(v.l).43^,jf5(v.2)W3^.5f^(v.3)  ’ 


(1-51) 


For  the  ^q's  (n  »  l,2,.,.f,2v)  we  then  use  (1*50) > 
For  n>2vel  equation  (1*47)  yields: 
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^  "^^n-l  ^1  ^  )  +^^-3  ^3  ^  ^“‘^*‘3 )  ■*‘^a^4^4  ^  n- v»4 )  = 

=  »a.2V-l*'6<»*''>*‘fi-tv./t(“‘''*l>»‘n-Svr.J*’s(“*''*^>  •  (^*52) 

Remembering  that  1^  that  this  ec^uation 

le  valid  for  all  Integral  values  of  n^  including  (1-4S)  as  a 
speeial  casSi  The  supplementary  initial  conditions  are: 

b^j=  b^g=  0  i  bg  as  in  (1-51)  ^  (1-53) 

Rutting  n=o  in  (1-52) »  we  get:  b^fQ(0-v)+b^^f^(o-v»4)  =  o  j  since 
f^(0«v)  =  0,  according  to  (l«38)i  we  obtain  b_^=  0^  By  repeated 
applications  of  (i-52)  we  also  find:  b^g=  b_g=*i*i.=  o*  ThuSj 
conditions  (1-53)  are  enough  to  satisfy  all  requirements  for  the 
coefficients  b^  of  S(x). 

The  process  so  far  has  left  undetermined «  This  does 

not  cause  any  difficulty.  Looking  at  (1-52)  observe  that  if 
it  can  be  satisfied  by  the  set: 

^0’  H*  ^2’  ^2v+l’  ^2v+2 . 

it  esui  also  be  satisfied  by  the  set: 

^o’  ^1'  ^2v+i*^®o-  ^2Vt2‘*’^H*  ***** 

where  k  is  a  constant*  The  reason  is  that  the  set: 

®«2v— 1^  0»  s^2v^  •».•»  Sq* 

satisfies  the  homogeneous  part  of  (1-52)*  Actually  for: 

n  =  2v*l«i  ,  lj„s  b2v*j+,=  (1-5*) 

equation  (1-52)  becomes: 

\^o  ^  ^  ^-2^2  ^  *®m-3^3  +B^.4f  4  (a*v-3 )  - 

*  S^o^®'''^**'^^‘*‘*m-l*’l^®*^^''’®B-2^2^®''’'^^^^  •  (1^55) 
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Its  homogeneous  part  is  satlsfiecl  by  the  set  a^>  ,  as 

a  comparison  of  (1*38) *(1*42)  *  for  n^+2v+l  (or  n*vsm+v+l)  ,  with 
(1*1) *(1*5)*  for  n+o’^+v+1,  reveals;  equations  (1*1) *(1*5)  with 
n+q=3n+v+l  give  the  coefficients  of  the  reeurrense  relation  (1*8) 
for  the  eoeffieients  a^  of  Ru^(x). 

We  can  choose  bgy^^=  arbitrarily.  The  simplest  choice  is 

“awl  =  So  *  0  • 

in  the  process  of  using  (1*52)  for  the  b^'s.  A  non*vanishing 
value  for  bg^^l^  means  simply*  that  S(x)*  and  consequently 
Rg(x) *  contain  the  additive  solution  b^^^^  R^(x),  which  can  be 
discarded  in  the  definition  of  the  second  independent  solution 
Rg(x)i  To  be  definite  then*  we  define  Rg^^)  ns  followst 

-  9^ 

R2(x)  =  Inx  R^(x)+x*''Z  b^#  j  |x|<fflin( |a|  *  |bj  )  ; 

n=so 

^2V+1®  *  ^0 

for  the  rest  of  the  b^’s  we  use  (1*52) * (1*55) •  This  definition 
leaves  no  affibiguity  as  to  which  particular  RgCx)  we  consider. 

jUialytic  Continuation  of  R^(x)  in  the  Right*Half  x*Planet 
We  use  the  Same  transformation  (1*15) »  as  for  R^(x).  In  analogy 
wltH  (1*30)  we  write: 

R2(t)  -  In  R^(t)+S(t)  *  s(t)  =  Z 

n=o 

For  a  moment*  designating  the  coefficients  of  (1*16)  by  T^(t), 
T^(t)>  we  can  write  this  equation  as  follows: 

T^(t)R"(t)+Tj(t)Rtt)+T2(t)R(t)  =  0  .  (1*59) 

Substitute  (1*58)  in  the  above  equations 


+  (t)dn(ff|)%+s3  =  0  . 

t(l-t)  2  _  1-t  1 

Since  Rj^(t)  itself  satisfies  (1-59),  we  obtain: 

T  ( t) 

tJt)a'^(t)4f.(t)8"(t)*Tg(t)S(t)  -  tut)'  w 

®  ^  ^  t^(l-t)^  ®  t(l-t) 


2T^(t) 

t(l-t) 


%'(t)  , 


Writing  f^Ct),  T^(t),  fgCt)  in  full  as  in  (1-16)  and  using: 

S(t)  =  i  ,  8,(t)  =  »,t“] 

n=G  »  1  msl  » 

we  obtain: 


2  t“*Vf(t.n.v)s(aa)''*^f  a„t®*’'*H(l-at)(t+b)(t.l)^(t*l)* 

n=0  0*0  * 

=  (aa)’'**!.  r(t,B+»+l)  ,  (1-60) 

m=o 

where  s 

F(t,in+y+i)  =  (iii+v+l)[2t^4(ah-4)t^-4ht^+4A(4h-2)t-2h]+ 

3 

♦[  (2^h)t\4(h*l)t^*(2-4h)t^+h]  =  2  F-(n+v+i)t‘^  (1-61) 

q=o  ^ 

f(t,n-y)  =  ^fa(n-v)t^  .  (1^62) 

q=0  ' 


The  funetlmis  appearing  in  the  last  equation  are  exactly  the  ones 
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that  appeafed  prevldusly  in  Rj^(t);  thus,  fg(x)  to  fg(x)  are 
defined  hy  equations  (1-18)  to  (1-24) *  On  the  other  hand,  the 
following  new  functions  were  introduced  in  (1-61): 


Fg(x)  =  -h(2x«l)  (1-63) 

F|(x)  =  (4h“2)x  (1-64) 

_  F2(^)  =  4x‘*^2-4h - — - - — {4-65) 

Fj(x)  =  -4hx+4(h-l)  (1-66) 

F^(x)  =  (h-a)(2x-l)  (1-67) 

Fg(x)  =  2x  i  (1-68) 

Aa  with  Rg(x),  we  call: 


®n"  Ho(2a)'^*^  g^  ;  g^s  1  ;  g^j^-O  ,  m  =  1*2,***  (1-69) 

and  write  (1-60)  in  the  following  foraj: 

HqI.  t^s_f(t,n-v)  =  F(t»m+v+l)  .  (1-70) 

°nso  ^  lasO  “ 

For  l^n^Sv  equation  yields  the  following  recurrence 
formula: 

6 

^  *'  ^  ,m  =  l,2,.,  (1-71) 

With  the  use  of  which  g^*  $2^  ••••»  g2y  nre  evaluated*  For  n=l 
it  yields: 

h  ^  I  -  .  (1*.72) 

For  n==?2v+i  equation  (l»70)  yields: 

®o^  ®2v^.l^0  ^  ^  *^2vh  ^  *®2v-1^2  f  »  *82v,p3^6  ^  ^  ^  ^  ^ 


1-16 

=  -h(2v+l) 

Since  =  di  aeccrdins  to  (1-18)  >  this  equation  aenves 

dimply  to  define  H^.  i^e* 

h(2v+l) 

Hvh  ^  ■^®2v-.2^3  ^  ''“2)  +g2^^3f^(  Vi.3)+ 


*82v-4^5  *®av-5^6  ^  ^ 


tl-73D 


For  n>2v4‘i  equation  (1-70)  yields  the  following  recurrence 
formula  for  the  remaining  coefficients  g^: 


^  X«n.2v-l-q 

qso  ^  ^  o  q=o  - 


(1-74) 


As  before,  the  procesa  has  left  s^yf^x  tnidetermlned*  it  will  be 
chosen  in  such  a  way  as  to  identify  the  present  solution  ^2^^^ 
with  the  solution  defined  explicitly  in  (1*57) •  Comparison 
of  (1-58)  With  (1-30) , (I-3I)  cuid  (1-15)  expresses  this  require¬ 
ment  as  follows: 

=  x'^'b^  Z  d  x’^  ,  (1-75) 

-  nso 

^o*^n  Substituting 


S(t)  =  Hg(2a) 


v+1. -V 


n=o 


where  as  in  (1-50)  we  can  write 


t  = 


x+2a  2a 


=  2L,(i+  2S_) 
—  2C' 


-1 


we  obtein* 


QP 


11=0 


(2a) 


2V+1  ^ 

n^o 


On*” 


(1-76) 


For  sufficiently  small  (x|  the  binomial  theorwn  yields: 


(!♦  ^  l+(v-n)|-  e 


(v-a) ( v»n*l) » . . . (v-n«m4l) 
m! 


'2a' 


(1-77) 


If  v«n  is  equal  to  a  positive  iateief>  it  reduces  to  a  finite 
polynomials  Assuming  v  itself  to  be  an  integer  (if  v^+l/2  the 
proGedure  remains  the  same) *  we  write  (1“76)  as  follows: 


H, 


fl+v  ^ 


v(v«l)  . 

»  . 


2! 


s«i(fe)Cis(v-i)  ^  * . ♦(fl)'"'!* 


■•‘®v+l'2a^  ^ 


v+2f 


‘*2  If  *  + 


(v+1) (v+2) 


/X  »2v+lr-s  /  ...%  X  .  ' ’L*'-!'""'' /X  s2.  T. 

*^2v+l'2a'  Ll-(v+l;  ^  +,*.,J+ 


aa 


a! 


■aa' 


2v+l 


^x+. . .  .♦Ox^'^^^+dgy^gX 


2v+2 


'+....] 


We  Immediately  deduce  that: 

H  -_!2 _ 

and  frcsD  the  coefficient  of  x  : 


(1-78) 


1-18 


+  ( *1 )  ^ 

v(v-l)  6 . , .  (v-i 


V4'3  v(v-l) 


4»  . 


(m-1)! 


2V 


(1«?9) 


Fof  v^+1/2  this  fortnula  is  modified  by  inclusion  of  some  addi“ 
tional  terms i  With  this  definition  (i«74)  can  be  used  to  yield 

-ftll-^ia  f on  n>  The  second  aolution»  defined  preeiaely - 

in  (1^5?) »  can  also  be  written? 

,v+l 


R2(x)=R2(t)=lnx  R^(x)+ 


H. 


y  *“[1+^  g^h  lt|<min(l,  |h|  )  (1-80) 
t  n=l 


and  so,  in  connection  with  (1-29) ,  can  be  analytically  continued 
Into  the  right-half  x«plane  shovna  in  figure  (1-1)  ^ 

Equation  (1-78)  can  be  checked  easily  for  low  order  functions 
(for  example  v=i,2j3)i  using  directly  the  definitions  (1-51) j 
(1-73) •  For  large  v  the  numerical  computations  of  Chapter  3* 

Pi^f  I,  checked  all  relations  and  proyed  them  correct*  in  Cases 
I  and  II,  both  series  in  x  and  t  were  used  to  evaluate  R^(x), 

the  results  were  identical.  It  was  also  found  that,  while 
the  series  in  x  for  both  Rj^(x)  and  R2(x)  had  the  same  rate  of 
convergence,  the  series  for  R^(t)  was  better  and  could  be  used 
farther  than  the  series  for  R2(t).  The  coefficients  g^^  Increase 
faster  than  the  Oj^'s.  As  the  order  v  of  the  functions  Increases, 
both  e^  and  g^^  increase  faster,  rendering  the  rate  of  convergence 
of  both  series  in  t  poorer*  In  any  case.,  the  series  in  t  could 
be  used  for  larger  |x|  than  the  series  in  x  and  were  able  to 
carry  the  computations  into  the  region  of  validity  of  the 
asymptotic  series  (especially  for  x  close  to  the  real  axis). 


ASYMPTOfiC  SOLUTIONS  R^(x)i  R^(x)  AROUND  X  =eo 


As  x^^^  tile  eoeffleletitfi  of  R^(x)  r'^(x)  in  equation 

(I  1*50)  vary  ae: 


_ 0 _ 

(x+a) (X4b) 


^  0(x*^] 


v(v+l) 


Point  x=6e  la  an  Irregular  singularity  of  rank  1  (10  ppi  58*77). 
We  follow  the  method  of  Chapter  III  in  referenee  10.  First*  we 
eliminate  R'(x)  hy  putting: 


y(x)v(x) 


(1*81) 


Points  x=-a  and  x=*b  are  branch  points  of  the  general 
solution  of  equation  (I  1-50)  and  in  the  cut  x*plane  with  branch 
lines  from  x=*a  and  x=*b  to  ©o  (along  the  negative  real  axis 
preferably),  the  functions  R(x)*  y(x)m^(x*$.)/(x*b) ,  v(x)  In 
(1*81)  are  analytic  and  single-valued*  We  obtain  for  v(x): 

v*'(x)+C2(y/y)+  ■**  W. r*»^iv/(x)->.{(y/y)^  e(y/y)^e(y/y)7  vV  "  ^  * 

(x+a)(x+b)  (xea) (xeb) 

c  v(v*l) 

*1*  ^  *  *y-*3v(x)  =  0  . 

But  y^y  =  »c/C2(x*a) (x+b) }  .  So,  we  get: 

y''(x)  +  q(x)T(x)  =  0  ,  (1*82) 


where : 


q(x) 


5/4  1/2 

-  e  s-*-*****!**^ 

(x+a)^  (x+a)(x+b) 


1+ 


+  f(i*  |)-i- 


-ijCl*  |)’®. 


V 


For  |x(>  max(  |a| ,  ffej )  ire  can  uae  the  binomial  theorem  to  obtain 
the  expansion: 


n=o 


A  few  manipulations  yield: 


^1^  ®  (1«I6) 


qgF  -V(v^l)-Cb 

^3"  (1-88) 


') 3 »  n^3»4»»»*  •(1*89) 


We  no>r  try  to  satisfy  (1*82)  with  the  fomal  series: 

v(x)=  e  0  ,  Ix|>max(  |a|  *  |b|  )  .  (1*90) 

Assximlng  h_gj=  0  (m=l,2,3» . . . . )  and  substituting  (1*90)  and 
(1*84)  into  (1*82)  we  obtain,  after  eancelllng  the  factor  e®* 

(10  pp.  88*77) : 


Putting  the  coefficient  of  x"^"^  equal  to  zero  wo  Obtain: 


1 


o  it 

in=o 


1-21 

(1-91) 


fox'  n^dili2,iii.  (for  n=-l,-2«>«««  it  is  automatieally  satisfied, 


since  q  h  0,  m=l,2i*iiO* 

-“in  “in 

2  2 

For  n=0  we  get:  #  0  or  a  +1 

For  n=l  and  since  a  +l  =  0,  we  obtain: 


Oi  That  is  eo^=  i,  dB2=-i. 


-2a>p+q^=  0  or  p  =  c/2®  ,  i.e.  p^=  -ic/2  t  pg=  ic/2  . 

From  (1-91) ,  replacing  n  by  n+1  and  using  =  0,  p  “  c/2®,  we 
also  obtain  the  following  recurrence  formula  for  the 

n+1 

2®nh^=  (n+e/2®)  (n+G/2®-l)hj^^2‘*‘  ^  ^^♦1-m  *  b^l*2M*»‘  •  (1-92) 

8=2 

Taking  h^=  1  we  can  find: 


( e/2® ) ( e/2®^l ) - V ( v+1 ) -eb 


and  So  on.  We  then  obtain  two  particular  normal  solutions  around 
x=oo  : 

R^(x)rv  JIH  e®^x“*^[l+2.  »  |x|  >  max(  |a|  ,  |bl  )  ,  (1-93) 

4 

where  R^(x)  is  given  for  1,  p^=  -lc/2  and  H^(^)  for  ®gS  -1, 
p^=  lc/2.  The  series  appearing  in  the  above  formal  representations 
are  normal  asymptotic  series  In  the  precise  sense  of  Foincare's 
definition  (9  pp.  168-174  444-443,  10  pp.  69-72),  In  the  region 
|x|>max(  |a|  ,  ib|  ),  Expanding  according  to  the  binomial* 

J^=  (1+  |)^/^(1+  =  (l+a/2x+...)(l-b/2x+....)  = 

=  l^c/2x<f . ...  ,  |x|>max(  )a| ,  lb| ) 
we  obtain  the  alternative  expressions t 
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R3  (x)  ^ 


(e/2«) (e/2»+l>-v(v+l)-cfe  * 

♦  ^  ♦•••][  1+c/ 2X*f  ••••]]■ 

2«d  * 


=  e®  V 


( c/2C8)  (  0/2®+!  )  -V ( v+1 )  «cfe 


>i*S 


2® 


63 

+  ^  +i 


]  . 


(1«94) 


which  will  serve  for  comperisofi  later^  Another  useful  remark  is 
the  following t  in  the  real  ease  (i.e*  a^h  real)  we  have  p^=  #2= 

=  aie/2f  where  the  bar  sipiifies  the  complex  conjugate^  while 
always  ®^=  5^^  i.  Then,  according  to  (l-85)-(l-S9)  i  all  Qji'b  are 
real  and  the  recurrence  formula  (1*92)  yielda  as  coefficients  h^ 
for  R^(x)  the  complex  conjugates  of  tjie  coefficients  of  R^(x)^ 
The  same  is  obviously  true  for  the  coeffiolents  in  (1-94) • 
This  mesuis  : 


R^(x)  =  R^(x)  ,  a,b  real  , 

(1^95) 

or  for  real  x: 

Rj(x)  =  R^(x)  . 

(1-96) 

The  critical  line,  or  Stokes  line.  Is  given  by  (10  p*  72) i 
Re(®x)  =  0*  If  x>x^^ix^,  with  ®  ^  we  havet 

ReC+l(x^+iXjL)  3  ^  ix^  =0  or  x^^  s  0  , 

So,  the  real  X-axis  is  the  critical  or  Stokes  line*  Any  solution 
R(x)  of  equation  (I  1-30)  can  be  expressed  asymptotically  aS  a 
linear  combination  of  the  above  formal  solutlpns  R^,  i.e* 

R(x)^r^  A^  H3(x)+A^  R^Cx)  (1-97) 

and  this  expression  holds  vmiformly  in  each  of  the  upper  and 
lower  half  x-planes  separated  by  the  Stokes  line  x^=  0*  The 
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eoefflicients  mayi  dt  course  *  change  values  from  One 

plane  to  another  for  the  same  solution  R(x)  (10  p.  73 )•  A  similar 
eomhination  expresses  r(x)  along  the  real  axis^ 

The  main  proiblem  of  this  research  is  to  find  explicitly 
expressions  lllce  (1-97)  for  the  particular  solutions  R|(x)i  R2(x) 
defined  in  the  preceding  section,  in  both  halves  of  the  x-plane 
and  along  the  real  axis.  This  is  the  subject  of  Chapter  2.  One 
last  remark  concerning  the  evaluation  of  Rj(x)  and  R^(x),  for 
complex  a,  b,  x  it  was  observed,  that  the  asymptotic  series  for 
R^(x)  could  be  used  earlier,  for  smaller  |xj ,  than  the  asymptotic 
series  for  R^(x),  The  order  v  did  not  affect  much  the  value  )x| , 
after  which  the  asymptotic  series  could  be  used. 


Fig.  1-1  Comformgl  mapping  of  on  to 

plane  according  to  (1-13) 


CMPflR  2 
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ASYMPfOtiC  EXFAKSIdNS  ^  R^(z)»  Rg(z)  FOR  LAROl  Izl 


INTRODUCTION 


We  are  going  to  employ  the  method  of  li  Ford  as  developed 
and  applied  in  Chapter  VIII  of  his  hooh^  reference  ?•  Starting 
with  r^(2),  Rg(z)  as  diefined  solely  hyj 

Z  a^z^]  »  \z\  <  ffiln(  |a|  *  |h|  )  (1-10) 

^  n=l  ^ 

R«(z)az’'^[l+%  h  z’^]  ,  |:Zl<min(  ia|  i  lh|  )  *  2v*l^lnteger  (l»ll) 

^  n=l  " 

'*‘^-1^1  ♦a^^jf  j  (n‘*‘0«3)  + 

♦a^^^f 2^(n+o’— A)  =  0  5  a_j^^  0  >  m  ^  1,2»3»»»»*  (1—8) 

around  z=0>  that  is  without  reference  to  the  fact  that  they  are 
solutions  of: 

R^'(a)+  r'(z)^^[1+  ^  •  ~^^]R(z)  =  0  ,  (I  1-50) 

(z+a)(z+b)  ^  ■  Z' 

we  are  going  to  obtain  their  asymptotic  expansions.  More  preci¬ 
sely.  in  a  manner  independent  of  the  previous  results,  we  will 
arrive  for  R-|^(z),  and  when  2v4>l  is  not  an  integer,  for 
expressions  like: 

R(z)  A,  R3(z)  +  A4  R4(2)  (1-97) 

Z=*po  ' 

and  at  the  s^e  time  determine  explicitly  the  values  of  the 
constants  A^,  A^.  The  procedure  will  also  yield  the  same  fimctions 


2«2 


R^(z)  and  R^(z)  in  exact  agreement  with: 


R^(z)  e^^z" 

4 


^  (c/2»)(e/^l)*»v(v+i)^efe  * 

^  2»  ®  Z^ 


%  1 

+  ; 

z^ 


-  |i 
4 


P3  -  I 


(!«' 


In  a  suhsequent  section  Ford's  method  will  he  extended  and 
applied  to  yield  the  required  results  when  2v+l  is  an  integer* 
Ford  does  not  treat  this  case  in  his  hook  hut  refers  to  it  as  a 
suhject  for  further  research*  Since  odd  integral  values  of  v 
appear  in  the  solution  R^j^Cz)  used  in  region  (2)  ^  figure  (1«2)» 
PiyRf  li  outside  the  antenna,  we  must  find  an  expression  like 
(1-97)  for  the  logarithmic  solution  R2(z).  Then  we  will  he  ahle 
to  obtain  the  linear  combination  R2|^(z)sA^j^Rj^(z)‘fAj^2R2(z) , 
permitting  the  evaluation  of  R^^Cz)  for  small  |z|  . 

We  are  going  to  make  use  of  two  fundamental  theorems  proved 
in  chapters  1  and  vi  of  Ford's  hook  and  refer  to  them  as  Theorems 
I  and  VI,  respectively.  They  are  as  follows: 

Theorem  _i_:  "  If  the  coeffloient  g(n)  of  the  power  series 


oo 

f(z)  »  radius  of  convergence  >0, 

n=o^ 

may  he  considered  as  a  function  g(w)  of  the  complex  variable  w  s 
=  x^ily  and  as  such  satisfies  the  two  following  conditions  when 
considered  throughout  any  arbitrary  right  half  plane  x^x^: 

(a)  is  single-valued  and  analytic 

(h)  Is  such  that  for  all  |y|  sufficiently  large  one  may  write 

|g(x+ly)|  <  Ke  ,  where  e  is  an  arbitrary  small  positive  qu^tl- 
ty  given  in  advance  and  where  K  depends  only  upon  x^  and  e,  then 
the  fiuxctlon  f(z)  defined  as  above,  Is  analytic  throughout  any 
sector  S  (vertex  at  origin)  of  the  z-plane  which  does  not  include 
the  positive  half  of  the  real  axis  and  f(z)  within  S  is  develo- 


pabld  asymptoticaly  as  follows i 


f (z)/V 


^  g(-n) 
n=i  z 


H 
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The  theorem  Is  supplemented  by  some  remarks  and  generaliza¬ 
tions;  we  shall  make  use  of  the  following: 

(I)  fhe  theorem  is  valid  if  g(w) *  besides  satisfying  condition 

(a)  j  is  sueh  that  we  may  write,  when  x>x,  and  )wl  is  large. _ 

A  ^ 

jg(w)|  ^  K  |w|  where  K  and  c  are  constants  of  which  the  latter 
may  be  positive,  negative^  or  zero^ 

(II)  If  f(z)  is  defined  by  a  power  series  of  the  type 

t(z)  *  2  g(fl)(2/l*)". 

n-0 

where  |i=0e  is  a  constant*  then  the  theorem  continues  as  before, 
provided  conditions  (a)  and  (b)  are  satisfied,  except  that  the 
excluded  ray  instead  of  argz  s  o  is  now  argz  = 


Theorem _yi:  "  Let  f(z)  be  a  function  of  the  eomplex  variable 
z  defined  by  the  series: 

h(n)z^ 


f(2)  •=  2 


n^G 


in  which  p  is  any  constant  (real  or  Complex)  and  In  which  h(n) 
may  be  regarded  as  a  funetion  h(w)  of  the  complex  variable  w  = 

;=  x-^iy  and  as  such  satisfies  ^e  two  following  conditions: 
h(w) 

(^)  p-p.— is  a  single-valued,  analytic  function  of  w  throughout 

Rw+p) 

the  finite  w-plane, 

(b)  h(w)  is  such  that,  when  considered  for  values  of  w  of  large 
modulus  lying  in  the  right  half  plane  R(w)=x>x^,  where  x^  is 
some  assignable  nvmber,  it  may  be  expressed  in  the  form: 

c,  Cp  Q  +6(w,s) 

h(w)=C„+  - 

0  w+p  (w+p) (wepe-l)  (wtp)(w‘^pvl),,,  (w4>p*t*8-l) 


g-4 

in  which  the  are  eonetants  and  lim  6(w,8)  =  Ot  8=d,i»2i>*« 
^  Iwl***© 

Then  I  for  valuee  of  z  of  large  modidust  the  funetion  f(z)  hai  the 
following  aepiptotie  expanaionet 


^  h^*nj 

f(z)A^p  X  ;  ‘«if/a<argz43ir/2  , 

n=l  r(p-n)  z“ 

f(z)Ay  e®z^"P  2=  ^  5  «fr/2<argz<ff/2  , 

n=  o  sr 

in  which  latter  developient  it  is  underetood  that*  if  z  s  pe^^* 
we  take: 

gl-P  ^  eC(l-p)(lhP*M)3  .  * 


Moreover,  if  in  (h)  the  quantity  may  be  regarded  ae  an  arbitra¬ 
ry  large  negative  number,  the  fxmction  f(z)  ie  developable 
aeymptotically  when  argz  »  %  tr/2  in  the  fozin  of  the  aum  of  the 
eeriee  in  the  above  two  expreazions  under  the  eame  intez^retatim 

for  z^“^*" 

The  following  renark  ca^  be  added  t  For  ^/2^argZ4^tr/2  an 
equivalent  aeymptotic  expansion  is: 

-  S,  ♦•**^*^  2  ^  ■'  •’'/2<«-8*<n/a  . 

n?=l  I  (p-n;z  n=o  z 


The  added  series  asymptotically  contributes  nothing  to  the 
expansion*  In  fact,  by  factoring  out  e^  we  havei 


f  (z)  rw  e®C«e"*  ^ 


^ 

n=l  i  (p-n)z“  n=o 


Now  for  *n/2<«urgz<ir/2  e”"  has  the  following  asymptotic 


expansion: 


e“*r*/  Oe 


2 

z 


z^ 


Then  the  aaae  Is  true  for 
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^6  li  ( 

Q*z  ^  Aj  0+  $  ♦  %  .  ,  “flp/a  <afgz<t/2  . 

n=l  r(p«n)z’^  ^  z^ 

Therefore,  the  two  expahslohs  are  aeymptotically  the  saae. 

SOLUTION  OF  THE  ASSOCIATED  DlFFlRENCl  EQUATION 

What  essentially  defines  Rj^(z)  and  Rg(z),  as  given  hy  (1«10), 
(l“ll) ,  is  the  recurrence  formula  (l“S)j  together  with  (l*l)*(l-5)* 
for  the  Coefficients  a^^  Following  Ford  (7  Chapt*  Vlll) >  we  repla¬ 
ce  the  index  n  by  the  continuous  variable  x  and  a^  by  the  general 
function  u(x).  That  l8»  for  x  =  n,  n  being  an  integer »  we  have: 

u(n)  s  a^  i  (2-1) 

Then,  the  recurrence  formtiia  (1-8)  tranefositB  into  the  difference 
equation: 

,4 

X  fg.(x+of»m)u(x-m)  m  0  • 

mso  ® 

Advancing  x  to  x+4,  i.e.  writing  x+4  for  x,  we  get: 


y  f^(x4>o*4-in)u(x+4-m)  =  0  .  (2-2) 

ni 

mso 

We  next  substitute: 

x+o  =  y  (2*3) 

u(x)  s  v(x+cr)  s  v(y)  (2*4) 

and  the  equatlpn  becomes: 

P4(y)^(y*’‘4)+P5(y)v(y*3)+P2(y)v(y+2)+p^(y)v(y+i)+ 

’•■pQ(y)v(y)  =  0  ,  (2*5) 

where,  with  the  use  of  (l-l)*(l-5)»  we  have: 


(y+3)“v(v-fl)  ]= 

sabC  (y^4)  (y+S)  ■^i(y+4)  -v( v+l)  ]  (2*6) 

P5  (y )  =f  3^  ( x+d^+3 )  =f  I  (y+3 ) = ( a+b )  (y+3)  ( y+2)  +c  ( y+3 )  *  ( a+b )  v(  v+l)  = 

=(a4b) (y+3) (y+4)*(a43b) (y+3)*(a+b) v{v+i)  (2*7) 

pg  (y )  ^2  ( x+or+2)  =f  2  )  -  (y+2)  ( y+i)  +a^*  V  ( v+i )  = 

=(y+2)  (y4»3)*2(y42)4a^*v(v+i)  (2*8) 

P|(y)=f3(x+Qr^l)=f3(y+l)s2a  (2*9) 

Po(y  )=f  4(  x+cr)  =f  ^(y )  =1  ( 2*10) 


the  last  expressions  for  the  p(y)'s  have  been  i^itteh  doim 
as  the  suitable  foras  fer  the  solution  of  eouatloh  (2*5)  by  the 
method  of  Laplace's  transformation  (7  chapt>  VIII >  17  Chapt^ 

XV  pp.  478*501,  18  Chapt.  Ill  pp.  57*|8>. 

In  order  to  be  able  to  apply  fheorems  I  and  VI  we  must  find 
that  particular  analytic  solution  v(y)  of  (2*5),  which  satisfies 
condition  (2*1),  or,  in  terms  of  y  and  v(y),  the  coadltlonst 

f(nvv^l);::aj^  for  all  integers  n  and  for  Rj^(z)  (2*11) 

T(n-v)=bj^  for  all  Integers  n,  2v4>l  not  an  Integer, 

and  for  R2(z)  •  (2*12) 

Equation  (2*3)  Is  a  linear  difference  equation  of  the  fourth 
order  with  polynomial  coefficients*  It  possesses  four  independent 
solutions*  For  the  general  theory  of  these  equations  and  espe* 
clally  for  certain  of  Its  results,  of  which  continuous  use  will 
be  made  In  this  analysis,  we  refer  to  Chapters  I,  III  of  refe* 
rence  18,  Chapters  XII  and  XV  of  reference  17  uid  Chapter  VIII  of 


2*7 

reference  7* 

The  method  of  Laplace's  transformation  Is  applied  by  assuming 

a  solution  in  the  form: 


(2*13) 


where  ^  is  a  line  of  integration  in  the  complex  t“plane,  suitably 


determined  later  and 


is  found  from  a  certain  differential 


equation  i  as  will-  be  seen  in-  the 
integration  by  parts  we  c^  get: 


(y+s)  Ut  =  I  tdt  =  [ty*®4r(t)3£- 

C 

-  [^ty*®^'(t)dt  (2*14) 


(y+s)  (y+s+1)  j  ty‘*‘®“^Mt  =  (y+8+1)  [  ty*®t ( t)  * 

,[ty*‘®*H'(t)]^H*  j^ty*®'*‘H"(t)dt  »  (2*is) 

Substituting  (2*13)  into  (2*5)  and  using  (2*14) , (2*15)  and 
the  last  expressions  for  the  p(y)'s  In  (2*6)*(2*10) »  we  easily 
obtain: 

Jjty*^Ct2^2(t)e"(t)*tf^(t)r'(t)+fj,(t)'Kt)3dt+[i(e,t)]t  =0  , 
where : 

^2(t)=abt^+(a+b) t^+t^=  abt^(t+l/a) ( t^l/b)  (2*16) 

fj^(t)=^2abt^*(a+3b)t^*2t^=  »2t^Cabt^+(a+b)  t+l]+ct^  (2*17) 

f^(  t )  =:*abv  ( v+1 )  t^*  ( a+b )  V  ( v+1 )  t^+[  a^*v  ( v+1 )  ]  t^+2a  t+l  = 

=  *v(v+l)abt^(t+l/a)(t-fl/b)+a^(t+l/a)^  .  (2*18) 

It  can  now  be  seen  how  tbe  last  expressions  for  the  p(y)'8 
in  (2*6) *(2*10)  are  suitable  for  the  formation  of  the  ^(t)'s* 


We  also  have  for  (1?  pp.  478*561 »  18  pp*  57*88) i 
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=  t(t) tyC(^|(t)  +  (y+l)fg(t)+t^|(t)]«t|r^(t)#*^^g(t)  .  (2*19) 

Since  (j)|(t)+tf|(t)  =  *2abt^(t+l/a)  (t+lA)‘*‘ct^+4abt^+3(a+b) t^+2t^= 
=  2abt^(t+l/b)  i  we  also  have: 

I(t,t)  =  abty‘*‘^(t+lA)^’Kt)[(y+l)(t+l/a)+2t]*t'(t)t(t+l/a)j  (2-26) 

We  conclude  that  (2*13)  provides  a  solution  of  (2*5)  If  'lf(t)  is  a 
solution  of  the  differential  equation: 

abt^(t4l/a)(t+l/b)t^t)*t(j)3^(t)i^'(t)+^^(t)t(t)  =6  (2*21) 

and  the  path  of  integration  t  is  chosen  so  that  I(t,t)  has  the 
same  value  at  both  extremities  of  the  path  when  it  is  open,  or  so 
that  I (tit)  returns  to  the  same  value  if  I  is  closed  and  t 
returns  to  the  se^e  point  after  describing  it. 

We  now  looK  for  the  behaviour  of  the  solutions  of  (2*21)  at 
the  vicinity  of  its  singularities.  The  equation  has  three  regular 
singularities  at  tj^=-l/a»  tgS-l/b,  t=po  and  an  irregular  at  t=0. 
For  its  solutions  around  t=0  we  put: 

♦  =  ,  g,.  1. 

n=o 


Then: 

§1  _ 
dt  ^ 

,bA2 

n^o 

n=o 

ifi  _ 

dt^ 

,m/tg 

h^o 

(n*^)  (n+g*i)g^#*^"^*2m< 

(n+p)g„t“*P’5  ♦ 

n=o  ^ 

♦  nV/* 

n=o  n=o 

Substituting  in  (2*^21)  and  after  the  factor  e^'^^t^  la  eanceiled 
we  obtain: 

Cabt^-t-Ca+b)  t+lH  %  (n^fp)  (n+p«l)g^t*^’*‘^*2m  2  (n+P)g  2  ^t“+ 

n=o  ^  n=o  n=o  “ 

42nil  g^t*^‘'^W2abA(a43b)t^2]£i  (n4p)g^t*'‘"2«a £ 
n=o^  n=o  n=o 

+{-abv(v+i)t^»(a+b)v(v+l)t^+{a^“v(v‘fl) )  t^4»2at+l]  %.  g^t*  =  0  . 

n=o 

Equating  the  eoefflclenta  of  n=0,l,2,3»‘ •  •  •  to  zero  we  have: 

For  t°  :  gQ^m^g^=  0  *  Or  m^+1  =  0  »  Or  ffl|@  1  ,  mg=  -1  * 

For  t  :  “2fflpgg+m'^g^+(a+b)m'^gi^+amg^-^g^+2ag^+gj^=  0,  or  aapsc, 

p  =  c/2a  ,  l.e*  Pj^s  *ic/2,  Pg-  lc/2.  So»  we  obtain  two  normal 
aolutlons  around  t°0: 

<kj(t)^  e^/S‘^°/^(l+g3^t+g2t^+,***.)  (2-22) 

i|r|j(t)  e“^/H^®/^(l%g^t+g2t^+*.*..)  i  (2-23) 

the  series  In  parentheses  are  asymptotic  as  t-^-O  (9  pp.  168-174 
444-445,  10  pp.  69-72).  the  theory  also  assures  that  are 

twice  differentiable  and  that  the  asymptotic  expansions  for 
i|fj,  tjj  (even  for  l^jj)  can  be  obtained  by  formal  differentia¬ 
tion  of  (2-22)  and  (2-23),  respectively  (10  pp.  58-77).  Now 
considering  'lfj(t)  and  ’^ijCt)  in  connection  with  (2-20)  we  observe 
that  if  t-»0  along  the  real  axis  we  are  going  to  have: 

I(ilr_,t)  si(ik--,t)  =  0  as  long  as  Rey> -1+ | Image |/2  . 

1 1=0  1 1^ 

Since  any  solution  of  (2-21)  is  a  linear  combinatiem  of  and 
^II 

If  Rey>-l+|lmagc|/2  ,  iC^Kt) ,  t]  o  (2-24) 

for  any  solution  ^(t)  of  (2-21)  and  as  long  as  t  goes  to  0  along 


g-lo 

the  real  t^axis.  This  statement  will  serve  te  fix  the  path  i  in 
(2-13)  so  that  to  provide  a  solution  of  the  difference  equation 
(2-5). 

The  iadiGial  equation  around  t^*-l/a  is: 
t(f).  (t) 

—  j.  ^ —  ] .  *  =  0  with  roots  #=0  and  2. 

afet^(t+l/b)  ^ 

t(^i(t)  ^ 

Around  t«=-l/b:  $ ( P-1 ) -[  — t  — r]  p  =  0,  roots:  P=0  and 

§  aht^(t^l/a) 


P(P-1)*E 


Since  tj^=-l/a  and  t2=“l/h  are  regular  singular  points  of 
(2“21)  there  exist  one  solution  'lfj^(t)  of  this  equation  around 
t^==l/a  and  another  '^^(t)  around  tg=-l/fe*  which*  in  the  neighbor¬ 
hood  of  the  corresponding  point  and  up  to  the  nearest  singularity 
can  be  expressed  by  the  following  convergent  series: 


4f^(t)  s  (t+l/a)^  [l+p2^(t+l/a)+p2(t+l/a)^+ . ]  (2-25) 

l^gCt)  s  l+q^(t+l/b)+q2(t+l/b)^+ .  .  (2-26) 

We  observe  that  ^^^(-l/a)  =  I'^Ci-l/a)  s  0  so  that  from  (2-20) : 
lCt2(t) it] 0.  Also  from  (2-20)  and  (2-26)  we  see  that: 

Combining  these  results  with  (2*24)  we 

conclude  that  there  exist  two  corresponding  solutions  of  (2-5)  in 
the  form: 

▼^(y)  s  PT  4^  <2'27) 

where  the  paths  and  Cg  are  shown  in  figure  (2-1).  The  dotted 
lines  represent  the  branch  lines  of  ’^^(t)  and  tg(t).  The  paths 
and  Cg  deformed  as  long  as  they  do  not  cross  branch 

lines  and  end  at  t=0  along  the  real  t-axls* 
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fhife  86  defined  solutions  v^(y)  and  v^^y)  ane  analytle  for 
(7  Ghapt.  VIII I  17  PP*  478-SOl.  18  Chap t.  Ill) : 

Rey>  *1+ |Image|/2  ^  (2-^29) 

It  will  he  showi  later  that  they  are  also  independent* 

In  order  to  obtain  the  third  and  fourth  independent  solutions 
of  (2-B)  we  put  (7  Ghapt.  Vlll) i 

v(y)  =  w(y)/r(y)  (2-30) 

and  substituting  in  (2-3)  we  obtain: 

.4 

X  Pjn(y)w(y+m)/r(y+m)  =  0  *  or 
m=o 

^4  (y )  w  ( y+4 )  ( y )  w  ( y +3 )  ‘^Pg  ( y )  w  ( y  *2)  *t‘P|  ( y )  w  ( y+1 ) + 

‘•‘PQ(y)w(y)  s  0  j  (2-31) 


where  with  the  use  of  (2-6)  to  (2-10) : 


f4(y)^P4(y)-ab{  (y^4)  (y*5)-2(y*4)-v(v^»l)  ]  (2-32) 

P3(y)=P3(y)  (y+3)=(a4b)  (y+3)  (y+4)  (y«^5)•(3a+5b)  (y+3)  (y+4)+ 

*  [  a+3P  -  ( a+b )  V  ( v+l )  ]  ( y+3 )  ( 2-33 ) 

?2(y)^2^^^  (y+3)-(y+2)  (y+3)  (y+4)  (y+3)-6(y+a)  (y+3)  (y+4)+ 

+[a^+6-v(v+l)](y+2) (y+3)  (2-34) 

yi(y)-P3^(y)  (y+i)  (y+a)  (y+5)^aa(y+i)  (y+2)  (y+3)  (2»35) 

-  (2-36) 


Following  similar  steps  as  before,  we  assume  a  solution  in  the 


form:  w(y) 


and  form  the  functions: 
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f^(t)  S  t^+1 

(2*37) 

^f{t)  m  (a«fh)t^«6t^+2at 

(2*38) 

fg(t)  =  *2aht^*(3a*«‘5h)  ACa^+6«v(v+l)]t^ 

(2*39) 

f-|,(t)  =  *2aht‘^+Ca»^3h*(a+b)v(v+l)  ]t^ 

(2*40) 

^^(t)  =  «ahv(v+l)t^  . 

(2*41) 

Then  t(t)  must  satisfy  the  fdllowing  differential  equations 


( t^+1)  t V''(  t )  -(j.3  ( t)  t )  t )  t%''  ( t )  ( t )  <r '  ( t)  + 

+^o(t)t(t)  =  0  . 

while  in  this  case  (1?  pp.  478*501,  18  Chapt.  Ill): 


(a*-42) 


3 


.  (2.*5) 


Now,  equation  (2*4t)  has  three  regular  singular  points  at 
tso,  tj-l,  t^s^i  and  an  Irregular  singularity  at  t=oo,  The  indl* 
Gial  equation  around  t=0,  according  to  (2*37) «( 2*42) ,  ist 
p(p-l)(p*2)(p*3)  =  0  with  roots  0,  1,  2,  3.  Four  independent 
regular  solutions  of  (2*42)  correspond  to  thes®  roots  and  ^y 
solution  of  it  can  he  expressed  as  a  linear  comhinatlon  of  these 
four.  These  solutions  are  of  the  general  fonn: 

’^g(t)^t  [fg(t)+fi(t)lnt+. . , ,  .4f^(t)  (int)®] ,  where  s=l,2,3,4  and 

^1^^^’  t..,.,  f^(t)  are  analytic  at  teo.  Then  t^^f  (t) 
vanishes  at  teO,  provided  that  Re(y+pg)>  0,  or  Rey>  0,»i,^2f *3° 

So,  if  Rey>0,  t^tgCt)  vanishes  at  t=5G  for  all  ♦gCt).  This  fact, 
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odmblned  with  (2»43)  means  that: 

If  Rey>  0  :  K’J'* |  t=o^  ®  (2*44) 

for  any  solution  4^(t)  Of  (2*43)*  The  result  Is  that  for  paths  t 
starting  and  ending  at  t=o>  we  will  have  [l(^>t)]j=  0  for  any 
Solution  ^f(t)  of  (2*43),  provided  that  Rey^O* 

The  indicial  equation  aroimd  t^=i  is: 

“$(1*1)  ($*2)  ( 3*3)  *^  t^f3( t)  ]/[  t^(  t«fi) 

with  roots  0,  li  2  and  “lc/2»  Around  tj^=  *1: 

Ip  ( p*l)  (p*2)  ( P*3 )  *  I  t^f^  (  t)  ]/[  t^(  t*i)  ^  ( p*l )  ( p-2)  =  0 

with  roots  d,  1*  2  and  p^@ic/2. 

As  before,  there  exist  two  solutions  ^^(t)  and 
equation  (2*42),  the  first  aroiaid  t=i,  the  second  around  ts«i, 
which  can  be  expressed  by  the  following  convergent  series: 

t3(t)=(t*i)^^‘^/^tl+h^(t*i)+h2(t*i)^^..,]  ,  |t-i|<l  (2*45) 

’^4(t)  =  (t+t)^®/^Cl+fl(t+i)+f2(t*i)^^,...]  »  |tH|<l  ,  (2*46) 

Then,  combining  with  (2*30) »  we  find  two  other  independent 
solutions  of  (2*5)  in  •^e  form: 

where  and  are  shown  in  figure  (2*2)  (7  Chapt.  VIII),  They 
start  and  end  at  tsO,  as  condition  (2*44)  requires,  enclosing  i 
and  *i  counterclockwise.  The  t*plane  is  properly  cut  by  the 
br^ch  lines  of  ^^(t)  and  ^/,.(t),  correspondingly,  shown  by  dotted 
lines  in  figure  (2*2),  The  so  defined  solutions  y^iy)  and  v^Cy) 
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of  (2*5)  are  analytic  for 

Rey>0  .  (2*49) 

We  have  found  four  Independent  solutions  Of  the  fourth* Order 
linear  difference  equation  (2*5)*  v^(y)  nud  Vgty)  defined  and 
analytic  for  Rey  >  *1+ |lmago| /2  and  given  hy  (2*2?)  and  (2*28)  *» 
v^(y)  and  v^(y)  defined  and  analytic  for  Rey>  0  and  given  hy 
(2*4?)  i  (2*48)  i  for  Rey *1+ 1  Image j /2  and  Rey<  0  we  define  the 
analytic  eontinuation  of  v^(y)i  v^(y)  *  respeeti* 

vely,  with  the  use  of  the  difference  equation  (2*5)  itself »  that 
is : 

Vg(y)  s  *[iP4(y)Vg(y+4)+p2(y)vg(y+3)+p2(y)Vg(y+2)+ 

+2aVg(y+l)]  ;  s=1,2,3j4  .  (2*50) 

We  see  that  the  so  defined  four  solutions  of  (2*5)  are  analytic 
in  the  whole  y*plane. 

We  are  now  going  to  obtain  expansions  of  these  functions  in 
terms  of  inverse  factorial  series  and  investigate  their  behaviour 
as  y**op  in  the  right  half  y*plane,  heretofore  indicated  as* 
y^«>  •  observe  that  V|(y)>  V2(y),  W5(y)=v^(y)r(y) , 

and  W4(y)sV4(y)r(y)  can  be  expressed  in  ^Uie  general  form: 

"s'y’  «  .  (8*51) 

"  •  s 

where  f_(t)  is  analytic  at  t=t„  but  not  at  t=0;  the  so  defined 
fveictions  of  y  are  analytic  for  Rey>y_,  (where  y_=*l*  llmagcl /2 
for  3=1,2  and  0  for  8=3>4),  Looking  at  equations  (2*21)  and  (2-42) 
we  see  that  t=0  is  a  regular  singular  point  for  f^(t)  and  f4(t) 
appearing  in  the  definition  (2*51)  for  w^(y)  and  W4(y) ,  but  an 
irregular  singularity  of  f]^(t)  and  f2(t)  associated  with  V||^(y) 
end  V2(y),  (Equation  (2*42)  has  also  an  irregular  singular  point 
at  t=QP»  but  all  the  paths  in  (2*51)  begin  and  end  at  t=Q 

O 
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without  approaching  the  vicinity  of  t^eo;  actually  *  this  is  why 
(2“5l)  defines  Ug(y)  only  in  the  right  half  plane  Rey>yg)*  The 
appearance  of  irregular  singularities  in  (3*21)  and  (3*^42)  is  due 
originally  to  the  fact  that  equation  (t  l«So)  for  R<j^(z)  and  Rgiz)* 
and  eonsequently  the  functions  themselves,  have  an  irregular 
singularity  at  z=66  (7  Chapt.  vili);  this,  in  turn^  leads  to  a 
non^nojmial  form  for  the  associated  difference  equation  (2*5)  (7 


Chapt*  villi  17  pp^  47S-501,  18  Chapt*  ill). 

When  t=o  is  a  regular  singular  point  of  fg(t)i  it  is  a 
familiar  fact  of  the  theory  of  difference  equations  that*  if  a 
sufficiently  large  positive  numher  m  is  selected #  then  Ug(y),  as 
defined  by  (3*51)  for  Rey>yg,  except  for  a  constant  factor 

0 

depending  only  on  co,  can  be  expressed  in  the  following 

o  S 

form: 


u,(y) 


»  r(y/«)  „(s) 

®  r(y/<D+$g+i) 


(2-52) 


where : 

^(s) 

n(y)  = 


1+Z 

npl 


(y+cDp .  ^cp)  (y^oP^+Sflo)  *  *  *  *  (y+ap_+ncD) 

S  S  P 


(2*53) 


is  an  Inverse  factorial  series  convergent  for  Rey>  y^  (7  Chapt. 
VIlI,  17  Pp.  485*487>  18  pp.  61*64).  We  can  apply  these  expansions 
Immediately  to  w^(y)  and  w^(y),  since  t=0  is  simply  a  regular 
singularity  for  f^(t)  and  f^(t).  Looking  at  figure  (2*2)  we  see 
that  and  can  be  deformed  into  the  circles  |t*ij  =1  and 
|t+i|=l,  respectively;  referring  to  references  17  (pp.  483*487) 
and  18  (pp.  61*64)  we  conclude  that  in  these  cases  we  can  take 
m  =  1,  v^(y)  and  v^(y)  are  defined  by  (2*43)  to  (2*48).  Comparing 
with  (2*31)  to  (2*53)  and  with  ©=1  we  obtain^  except  for  a 
constant  factor,  the  following  expansions: 

yj(y)  -  (i)yQ|y|/r(y+i*lc/2)  ;  Rey>  0  ;  i=e^^/^  (2*5^) 


▼4(7) 


(*i)y  n  ^yj/r(y+l+le/2>  ;  Rey>  0 


;  (2*55) 

^‘*  •  •  •  »  Roy  ^  0  (2^56) 


1+ 


i 


(4) 


y+i+ic/2  (y+l+lc/2)  (y‘«‘2+ie/2) 


*(2-57) 


The  inverse  factorial  Series 


converge  uniformly  fen 


Rey>  0.  Their  GOefficients  can  be  obtained  by  the 

method  explained  in  references  17  (Ghapt.  XV) ,  18  (Ghapt.  Ill) ^ 
or  by  direct  substitution  into  the  difference  equation  (2-5)  and 
application  of  the  method  of  undeteimiined  coefficients^  We  will 


find  them  explicitly  later  using  the  former  method*  For  the  ge¬ 
neral  theory  of  inverse  factorial  series  and  their  properties  we 
refer  to  references  17  (Ghapt*  X),  19  (Ghapt.  Vl  pp*  170*177)* 
The  constant  factor  for  v^(y)»  ^4(7)  was  selected  so  as  to 
give  them  the  form  indicated  in  (2-54)  to  (2-57) »  Which  will 
prove  of  convenience  later*  Hereafter  we  adopt  these  definitions 


for  v^(y)  and  V4(y),  instead  of  the  integral  forms  (2»47) » (2-48) . 
They  are  simply  constant  mialtiples  of  the  latter*  As  before, 


v^(y)  and  v’4(y)  are  analytic  for  Rey^O,  the  same  being  true  for 
^(5)  ^(4) 

and  iZ^yj  .  Their  analytic  continuation  for  Rey^Q  is  again 


provided  by  (2-50)  and  makes  the  so  defined  functions  v^(y)  and 
V4(y)  analytic  throughout  the  y-plane. 

When  teO  is  an  irregular  singular  point  of  f_(t)  In  (2-51) » 
convergent  expansions  of  the  form  (2-52)  and  (2-53)  can  not,  In 
general,  be  obtained.  However,  with  m^i  (2-52)  and  (2-53)  provide 
an  asymptotic  expansion  for  u_(y),  as  defined  in  (2-51) »  as  y«*pp 
in  the  sector  -Tr/2+e  <argy  <C''^/2+t  (e>o  and  arbitrarily  small) 

(7  pp.  309-318,  8  pp.  70-74,  17  pp.  457-459,  20).  This  statement 
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Is  based  upon  a  tbeorem  Gontalned  in  nefeifenGe  8  (pp.  70-74)  .  It 
Is  due  originally  to  Norlund  and  its  proof  Gan  be  found  in 
Ndrlund's  original  paper  in  "Acta  Mathematlca",  reference  20.  We 
give  the  theorem  in  full,  as  contained  in  reference  8  by  Ford, 
because  It  serves  to  clarify  many  points  regarding  the  solutions 
of  difference  equations  with  polynomial  coefficients,  which  we 
have  been  making  use  of. 


NQrlmdl.s_„Theoremi  "Ol van  the  linear  difference  equation 
k  ‘ 

2]  F.  {x)u(x.-l)  =  0  >  (I) 

i=o  ^ 

where  the  coefflGlents  are  factorial  series  of  the  form: 


Pl(x)=Co 


.(Di 


.(1) 


x+1  (x>fl)(x<f2) 


i  s  o*l,2,..».*,k 


all  of  which  converge  throughout  the  right  half  of  the  x^plane* 


Suppose  first  that  the  roots  a^,  a^* 
riitic  equation 


a^  of  the  character 


{k>„ 


°0^^^  0  ,  c^^^^  0  (III) 


are  distinct.  Then  there  exist  k  solutions  u^^,  Ug, . . 

(I)  such  that  throughout  the  sector  -tt/S+s  ^  argx 
arbitrarily  small  and  ^0)  we  havet 


r(x-fi) 

r(x-Pj+i) 


f 


of 


(IV) 


where  is  a  constant  and  i^jCx)  a  factorial  series  of  the  form 
Indicated  in  (II) .  In  case  (III)  has  multiple  roots  and  a^  If  an 
n^fold  root,  two  cases  are  distinguished: 


(1)  Sj  is  at  the  same  time  an  (n-p)-fold  root  of  the  equations 


k 

z  ° 
8=0 


(a)  k^s 


=  0  ; 


n»l  • 
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(2)  These  eoniltidhs  are  hot  fullfiled. 

In  (2)  no  asymptetiG  development  exists  of  the  form  (IV).  In  (1) 
there  exist  n  linearly  independent  solutions  Ug(x)  (ssl,2f * « .* »n) 
such  that  when  -fr/2^‘e<  argx<ff/2-e  we  have: 


I  a  a  •  a  t 


n 


Ug^a|$g(x)  ;  S=l,2,. 

where : 

r(x+l)  V  r(x-fl)  vh  r(x+l) 

®  °  r(x«p^4^i)  ^s  Rx-p^+i)  ^p2  r(x*p^»»-i) 


the  expressions  being  developments  of  the  form 

(11). 

If  some  of  the  roots  of  (III)  are  zero  or  infinite,  it  is 
necessary,  in  order  to  obtain  a  system  of  fundamental  solutions* 
to  use  a  series  of  substitutions  of  the  form: 


U(x)  s  [Rx)] 


'  (x) 


''(X) 


(u^) 

'"(x) 


and  dete«fline  so  that  the  difference  equation  in  w 


(Pp) 

(X) 


shall 


have  a  characteristic  equation  containing  at  least  one  root  which 
is  finite  and  different  from  o.  It  is  always  possible  to  determiine 

in  but  one  way,  a  series  of  numbers  . .  ?uch  that 

the  total  nvffliber  of  roots  which  are  finite  and  different  from 


zero  in  the  corresponding  characteristic  equations  thus  obtained, 
is  exactly  the  order  k  of  (I).  If*  whenever  a  multiple  root  occurs 
in  one  of  these  characteristic  equations,  the  corresponding 
conditions  under  (1)  are  satisfied,  then  there  exist  a  system  of 
fundamental  solutions  of  (I),  each  of  which  is  asymptotically 
represented  within  the  sector  »Tr/2+6  <  argx <  TT/2’^e  by  a  series  of 


Up  X 

the  form:  Sj  $g(x)  .  If  no  multiple  roots  occur  in  one  of 

these  characteristic  equations*  the  solutions  have  the  simpler 
form: 
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>r  X  nxH^l) 

r(x) 


Exceptioas  occur  when  some  Of  the  numbers  are  not  Integers* 

Then  the  coefficients  in  the  above-mentioned  difference  equations 

(U:  ) 

for  are  no  longer  developable  in  factorial  series  of  the 

form  (II).  Suppose  Up  Is  equal  to  a  rational  fraction  p/q*  We 
put  xspz,  w(x)=.v(z)  and  derive  from  (1)  a  difference  equation  fer^ 
v(z),  thus  demonstrating  the  existence  of  solutions  expressible 
asymptotically  in  the  forasj 

r(U)  "f"  ^*<*/*’> 

Here  we  are  dealing  with  the  difference  equation  (2-5)* 
Dividing  by  P^(y)  and  making  appropriate  change  of  variables, 
like  y^-fsxj  v(y)=u(x),  we  can  reduce  it  to  the  forms  (1)*  (II)* 

We  can  then  easily  verify  the  statement  made  on  page  2-16* 

Returning  now  to  (2-5l)  we  observe  that  f^(t)  and  f2(t)j 
appearing  in  the  definitions  of  Vj^(y),  V2(y)  for  Rey>' -If  |Imagc]/2, 
have  an  irregular  slngxilarity  at  t^o*  Then,  according  to  the 
preceding  discussion,  (2-52)  and  (2-53)  with  ®=1  provide  at  least 
asymptotic  expansions  for  v^(y)  and  V2(y)  when  •  Referring 

to  figure  (2-1)  and  to  the  definitions  (2-25)  to  (2-28),  we  obtain, 
except  for  a  constant  factor,  the  following  expansions: 

. **  r(y)  ^(1)  .  . 


Vj^(y)  =  (- 


r(y) 


y(y+i) (y+2) 


'"(y) 


Vy)-(-l/b)y^^Q(y,  .  (-iA)yD(y,/y 


(1) 
(y)  y 


y+3  (y’»'3)  (yf4) 


(2-58) 


(2-59) 


(2-60) 


^ . 


(2) 

t  s  ~ 

(y> 


(y+l)(y+2) 


2»20 

(2#61) 


All  these  expressions  are  valid  in  the  seotpr  ■^tr/2+e  <  argy  <ff/2»6. 

^(1)  .(2) 

The  inverse  faetorial  series  *  ^^(y)  least  asympto^- 

tic  as  in  the  indicated  sector^  Their  eoefficients 


F=li2i  can  he  Ohtalned  hy  both  methods  mentioned  previously  for 

only  difference  being  that  the  process  is  now  at 
least  formal.  As  we  shall  see,  it  will  not  be  necessary  to  find 


them  explicitly.  Hereafter  we  adopt  the  above  definitions  for 
v^(y) »  VgCy).  They  are  simply  constant  multiples  of  the  Integral 
forms  (2»27)  j  (2-^28)  of  the  analytic  functions  v^(y)»  v^Cy). 

The  linear  independence  of  the  so  obtained  solutions  V|(y), 
V2(y) >  v^(y)»  v^(y)  can  be  shown  by  referring  to  Norlund's 
theorem,  or  by  making  use  of  a  general  theorem  given  in  reference 
17  (p*  360),  based  upon  the  behaviour  of  the  solutions  v^(y)  (ss 

D 

-l»2»3i4)  as  y^j^  »  which  we  now  know. 

We  have  obtained  for  v^(y),  v^(y)  two  convergent  and  for 
Vj^(y)j  VgCy)  two  at  least  asymptotic  expansions  in  terms  of 
inverse  factorial  series  valid  in  a  half  plane  limited  to  the 


left.  It  is  a  fact  of  the  theory  of  linear  difference  equations 


that  when  |y|  is  large  (without  any  limitation  to  the  left) •  the 


corresponding  expansions 


(s 

(y) 


though  not  necessarily  convergent 


any  morei  are  in  any  case  developable  asymptotically  in  series 
of  the  sMue  fonn  (7  pp.  3Q9»3l8f  17  pp»  457-459).  Thus*  we  have: 


v^(y)  -  (-1/a) 5^ 


H’fi(y) 
y(y+l) (y+2) 


(2-62) 


Vg(y)  =  (-i/b)y 


i+e2(y) 

y 


(2-63) 
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Vjty)  =  <l)yA|y|/r(y*l-ie/S)  i  (2-64) 

V4(y)  *  (-DJ^nlyl/rCy+l+lc/a)  ;  .i±e*^«/2  ,  (2-65) 

where : 


lla  4f,(y)  =0  ,  8=1,2 


where  the  upper  sign  should  be  used  for  8=3 » 


(2-66) 

. . *  8=3,4  (2-67) 

the  lower  for  8=4i 


ASYMffOTIG  KPANSIONS  OF  R^(z),  R^(z)  FOR  U^OE  |z| 

We  are  n(Dw  In  a  position  to  apply  these  results  and  fheorems 
I  and  VI  to  obtain  the  asymptotic  expansions  of  R|(z)j  ^^(z)  (for 
the  latter  when  2v+l  is  not  equal  to  an  integer),  solui-ions  of 
the  differential  equation  (1  1-50).  As  was  stated  in  the  preGeding 
section,  we  are  looKing  for  the  particular  solution  v(y)  of  (2-5) » 
which  satisfies  the  conditions 

v(n+o)  =  for  all  Integers  n  (2-68) 

and  where  for  R2^(z)  we  use  o=(7^=v+l,  while  for  (for  non- 

Integral  values  of  2v+l)  osOgS-v,  We  have  obtained  four  indepen¬ 
dent  solutions  of  the  fourth  order  difference  equation  (2-5) •  Then: 

v(y)  =  Av^(y)+Bv2(y)’*-Cv^(y)+Dv^(y)  ,  (2-69) 

where  A,  B,  C,  D  are  cons  tan!;:,  it  is  now  obvious  that  with  cr=v'+i» 
or  cr=-v  (if  2v+i  is  not  an  integer)  and  the  following  fovu’  initial 
conditions: 


V(0+0')  =  v(0’)  =  1 

(2-70) 

v(«i+o)  S  0 

(2-71) 

V(»2>fO')  =  0 

(2-72) 

v(-3+o')  =  0  , 

(2-73) 

.vCy+o)  will  satlsi 

■  y  the  c-oadltlons  (2-68) 

for  all  Integral  values 

df  n.  The  proof  is  as  follows:  For  n=0,li 2,3 » • • • •  aad  by  virtue 
of  the  above  relatioas,  v(n+o')  will  satisfy  the  same  reourrence 
formula  (1-8),  whieh  the  eoefflclents  a^  satisfy,  sinoe  for  y^<fcr 
(nsO,! , 2, . . . )  the  dlffereaee  equation  (2-5)  reduces  to  the  re¬ 
currence  formula  (1-8),  ir  fact,  (2-5)  was  derived  from  (1-8)  ifi 
the  preceding  section*  Thus  we  have:  v(n+cr)sa^,  n=0»lj2,3,  *  *  * , 
verifying  (2-88)  for  such  n's.  Oh  the  other  hand,  for  ys-4+qf 
(2-5)  yields: 

P4  ( -4+b)  V  ( &)  ( -3+d')  V  ( -l^o’)  ( -2+cr)  y  ( -2#cr)  ^2av  ( -3+ d)  4^?  ( -4+or)  =0. 

Using  the  last  three  conditions  (2-71) -(2-73)  we  get: 
p^(-4+0')-fv(-44-cr)  =  0. 

According  to  (2-6):  p^(-4+cr)sab(cr-v-l)  (a^-v) .  iinCe  cr=v+l  or  crs*v, 
we  see  that  in  either  case  p^(-4+cr)=0  and,  consequently,  v(»4+cr)  = 
=0  (q'=c^  or  o^Og) .  Using  (2-5)  and  v(-^n*a')=0  for  nal,2,3»4  we  see 
also  that  v(-n+o')=P  for  n=5,6,7,--*  •  Therefore,  v(n+o)=0=a^  for 
n=-l,»2,»3,-4,-5, • • • • •  and  (2-68)  is  verified  for  all  integral 
Values  of  n,  as  required*  The  initial  conditions  (2-70) -(2"73) 
serve  at  the  same  time  to  determine  the  coefficinrs  A,  1,  C,  D 
in  (2-69)* 

As  a  consequence  of  all  these  results  we  can  write* 

R,(a)  =  z' X  ^  ^  ^(n+cr)?  ,  where  of=?/+l  is  used  for  R,(z) 

i  n=o  ®  n=o  ^ 


a-*  23 

anii  wiien  2v+l  is  not  an  Integefi  f=-*v  fon  Rg(z).  Using  (2“69) 
for  vCn+o)  we  obtain} 

R|(z)=  A(-z/a)^£  G^(n)  (-z/a)*^+B(^z/b)^2  ajn){^z/b)% 

2  ft-o  n=b  - 

^G(iz)*^2.  ^(n)(iz)^+l)(-iz)^Z  G.(n)  (*lz)^  *  (3-74) 


i  e*  Gg(n)  when  eonsidened  as  fynetions  of  w=x+iy  ane 
given  by: 


&|(w)^|(w+cr)/(-l/a)' 


Gp  ( w)  =Vg  ( w+o) /( -l/b ) 


(w+o)  (w+0+l)  (Mi*a*2) 


(3-75) 


(3-76) 


iim  e. 


=  lim 


lW|«*-e 


=  0  » 


G,  (w)=v^  (w+o)/(i)'''*’^s  Q  rL«:i/r(w+0+l-ic/3) 


(2-77) 


G^(w)=V^(w^qr)/(-i)'^*^8  Q  (w+0+l^ic/2) 


p,(3) 


A3) 


'\w|-*.eo  w+a+l-^ic/3  (w+0+l-ic/2)  (w+a+2-ic/2) 


(2-78) 


(2-79) 


(4)  g^^^  g^^^ 


the  last  expressions  in  (2-73) -(2-78) ,  as  well  as  (2-79)  and 
(2-80)  being  obtained  with  the  use  of  (2-62) -(2-67) ,  With  Re(w+0) 
yo,  or  Rew^ -a,  (2-79)  and  (2-80)  are  convergent  series.  The 
functions  Gg(w)  are  analytic  for  all  w  Just  as  the  functions 
Vg(w‘fa)  are.  Furthermore,  according  to  ( 2-78)- (2-80) ,  G3^(w)  and 
G2(w)  satisfy  all  the  conditions  of  Theorem  I  subject  to  the 


2-24 

remarks  i)  >  il)  *  while  G^(w)  aaid  CJ^(w)  satisfy  all  the  conditions 
of  fheorem  VI  with  P2=dfl»te/2  and  According  to  the 

ratio  test  and  with  the  help  of  (2-75) “(2-80)  we  see  that  the 
first  series  In  (2-74)  is  convergent  for  |z|  <  laj  ,  the  second  for 
ijzl  <  Ihl  ,  while  the  last  two  are  entire  functions*  The  first  two 
can  be  expanded  asymptotically  for  large  izl  by  applying  fheorem 
li  the  last  two  by  applying  fheorem  VI.  With 

z  =  Izje^^  «fr^(^^tf  (2-81) 

we  have  using  fheorem  1,  1) ^  11) j 

Al-z/a)*^!  G.  (n)  (-z/a)^^-A(-z/a)^i  G.  (-n)/(-z/a)^^ 
n=o  ^  n=l  ^ 

=-Az®X  »  for  ^  ^  arg(-a), 

nsl  ^ 


or  using  (2- 


A(-z/a)^2 

n=o  nsl  z 


^f>o  Avj^(-n+or) 


Similarly  with  the  use  of  (2-76) s 

Bv«(-n+cr) 

B(-*A)^Z  &-(n)  (.s/b)“/V  Vl 


n?=o  2  '  ^1  z^ 


,  -ft<  ^  ^TT  * 

^  fi  arg(-a)  .  (2- 


-Tr<(^^TT  , 

^  ^  arg(-b)  .  (2-83) 


For  the  third  term  in  (2-74)  we  apply  Theorem  VI  with  the 
accompanying  It  remark  and  with  p=o+l-lc/2.  Also  here 

l=e^^/^  and  so,  with  (2-81),  -Tr/2  ^arg(lz)^  5tt/2.  Thus! 


CG,(-n) 

C(lz)^i  G,(n)(lz)^/N; -(lz)®Z  »  TT/2<arg(iz) 

n«o  ^  ^ 


/V- 


n=i 


or  0<^^TT 


n=l  (Iz) 


1+: 


(5) 


n^l  (Iz)- 


r]  ; 


“fr/g  < afg(lz) 4  ff/s  oi*  *ft:<(|»4o  i 
AccoMing  13  the  theoz^em  in  the  second  expansion  we  must  take: 

Since  in  this  case  “Tr<^4  0i  we  must  take  ’^=<^^fr/2>  hecause  only 
in  this  way^  wa  can  have  -tt/2-^- agfeement  with 
Then  we  obtains 


Using  (2-77)  and  (2-8A)  we  finally  obtains 

fr^  «  ( ‘^n+o) 

C(iz;)°^2  U^(n)(iz)”A^  -Z^Z  — ,  0<(f)<fTi 
n=o  ^  n=l  ^ 


nsl 


n  Cv»(-n+a) 


ah 


z 

hil 


6. 


(5) 


n^ 


For  the  last  tei*m  we  haves  ^^ise"-  ' -  , 
4tr/2»  pso+l+ic/t.  fhens 


(2-84) 


-fT<.^  4  0(2-85) 
‘3Tt/2  4arg(-iz)4 


„  eo  DO,  (-n) 

D(.tz)®2  <}4(n)(-U)"<V.(.l8)‘'l  -4-— 

n=o  ^  nil  (-Iz)’^ 


DG I  ( -n)  j 

^  -(-iz)<^2  ^(,iz)  W^"(-iz)’‘^-’ 


nil  (-Iz) 


n 


-3Tr/2<  srg(»lz)<  -tt/2 
or  "^IT  <  <  0 


li*z 


n=l  (»iz) 


-w/2^arg(-lz) ^ tt/2  or  o4^4tt  . 

In  this  cases  (-iz)*^"^°^-??  g(»ff-lc/2)  (in|z|+l^)  ^  where 

»TT <t^arg(*lz)4  TT.  With  04^4^  must  take  t=^»Tr/2»  resulting 

in  <-Tr/2  4'l' 4^/2»  in  agreement  with  -Tr<lT4w  Thuss 

(„l2)^<y«lc/2_^-a-lc/2^(-e-lc/2)  (-ln/2)_  jj^®^Ac/2^^j^)-a^-Trc/4^ 


This  result  aiid  (2^78)  finally  yield: 


D(^lz)  02.(n)(“iz)  rvj 
n=o 


n=l  z 


n 


rj  »z 


fDe 


Ilfil  z 


n 


n=l  (» 

j  -ir  <  ^  <  0 


(4) 

““SJ  . 


Comfeinlng  (2*^82) ,  (2-83)  *  (2“85)  and  (2-86)*  substituting  In  (2*74) 
and  malEing  use  of  (2-69)  we  obtain: 


^  ^(-n+cr) 


Rt(z)a;«z  2 

«  n=i  z 


^De 


-iz  ^* 


(4) 


n_ 


.J»|  2W/2[i*2,  -^]  .  .«<♦<' 

n=l  z^  nPl  (iz)^ 


n=l  (-iz) 
(3) 


.  0<f<fr 


.^••1  /  j  _ 


(3) 


•3^ 


nsi  z 


+De 


-iz^- 


2  7^^ 

n*l  (- 


nsl  (Iz)’^' 

*i  ,  (j>=o  5  ^^arg(-a,*b)  (2*87) 


The  ray  Is  a  branch  line  for  ^^(z),  ^2(2)  and  is  exclu¬ 
ded  from  the  above  expansions.  For  large  |z|  we  can  draw  the 
branch  lines  1  starting  from  z=-a  and  z^-b»  along  ^^ir.  Also 
according  to  (2-68):  f(-n+cr)  =  a_j^  =  0,  n=lf2f3>*«*i  so  that  we 
finally  obtain; 


( z)  Pe" 


e*^^  z* 


rv  ce”’^®/^ 


e  ”  z 


00 

Ci*i  =ft-^3  .  0<K« 

Ml  (.Iz)® 

«  8<’> 

u-z,  J  .  •«<:<o 

Ml  (iz) 


2.27 

,  ,  »  .(?) 

■  ml  (Is)" 

M  »(*) 

♦£«'*"/*  -ia--r3  .  ♦  e  0  .  (3*88) 

n=l  (*l2)^ 

For  R^(z)  we  use  a=vH»l2  for  Rg(z)  (when  gv+l  is  not  an  Integer)  ^ 

#=-Vi 

In  order  to  complete  the  prohlem  we  are  going  next  to  detenni 
ne  explieitly  the  GoeffiGients  end  identify  the  above 

expansions  with  the  previously  determined  asymptotiG  solutions 
R^(z)j  R^(z).  Mtieipating  the  forthGOming  proof,  we  can  state 
that,  in  an  entirely  independent  manner,  we  have  arrived  at  the 
required  expansion  (1»97)  forR^(z),  R^Cz),  in  complete  agreement 
with  the  results  of  the  preceding  chapter  and,  at  the  same  time, 
have  Obtained  explicit  relations  for  the  precise  evaluation  of 
the  coefficients  of  the  linear  relations* 

For  g^^^  we  use  the  method  explained  in  references 

XI  n 

1?  (Ghapt.  XV)  and  18  (Chapt.  Ili) *  mentioned  on  page  2*16.  We 
first  deternine  the  coefficients  h^^  of  the  solution  '^j(t)  around 
t=l  of  equation  (2*42).  'k^(t)  is  given  by  (2*45).  We  put  t-l=z, 
t=z*l  and  substitute  in  ( 2*37)  ■^( 2*42)  s 

z(z^2l)  ^*4  *[  (a+b)  (z'»-l)'*6(z+l)+2a]i*t  ‘»-[ab(z+l)^*(3a+5b)  (z+l)  + 
dz^  dz^ 

2 

-fa^-*>6*v  ( v+l )  *[  *2ab  { z+1 )  ♦a*3b*  ( a+b )  v  ( v+l )  3|^  * 

dz^ 

*abv(v*l)t(z)  =  0  .  (2*89) 

Assume  ’^^(z)  «  h^Z"'*'^  and  substitute  in  (2*89).  After  multi* 

h*Q  ■  ^—6 

plying  it  by  z-  ^  we  obtain? 

£  (n+p)  (n+p*l)  (n+p-2)  (nep*3)h^z“(z+21)*  J  (n^P)  (n+P*l)  (n+p-2)hj^’ 
n=o  ”  n=Q 


‘  z’^l(a+b)  z^+2Cl(a+ib)  »3lz+e*4l|+  2  (*♦# )  (n+fwi)]i,z^ral3Z^'«(3a‘fSb“ 

*2iab)  z^+iaef^^vCvH*!)  «i(3af5fe)  ]z?+  ^  (nfg)h;^Z*^f2abz^+C  (a+b)v(v4‘l) 

J  nSo  “  ^ 

-a*3b+2iab]z®f- X  bi^z*^abv(v+l)z^=  ^  z^hjjf  (z*n+#)  =  0  » 

^  n=:d  h— 6 

wb&re : 


f  (z,n+p)  =  (ri+P)  (n+p^l)  (n+p-2)  (n+P-3)  (z+2i)*(ft+P)  (n+p“l)  (n‘»‘$»2)  ‘ 
‘^a+b)  z^+iCi  {a+b)*3]z‘fc«6i^  +  {n+p)  (n+p-l)^abz^»(3a‘»>5b“2iab)  z^+ 

‘f  t  ac+6^v  ( v+1 )  *1  ( 3a+Sb )  ]  zj  +  (n+$ )  ^2ab  z^+[  ( a+b )  v  ( v+1 )  *a'“3b+2iab3  z^j 
«abv(v+|)z^=  f^(n+p)+f^(n+p)z+f2(n+p)z^+fj(n+$)z^  * 

where  t 


f^(x)=:3c(x-l)  (x«2)  C21(x*3)-c46i]  (2-90) 

f|(x)=x(x-l)  ^(x-2)  (x-3)+2[3-l(e*b>  3  (x-2)*«'ac+6-v(v+l)- 

-i(3a+5b)^  (2-91) 

f2(x)=x[-(a+b)  (x-^1)  (x-2)-(3a+5b-2iab)  (x-l)  +  (a+l>)v(v+l)- 

-a-3b+2iab]  (2-92) 

f^(x)=ab[x(x-l)+2x-v(v4-l)]=ab(x-v)  (x+v+1)  .  (2-93) 

Reeurrence  relation: 

3 

J  fjjj(n+p-m)h^_jjj=  0  ;  hQ=  1  ;  h  .=  0  ,  J=l,2,,.,.  ,  (2-94) 

nf=o  "  - 

Indlclal  equation;  fQ(P)=PO-l)  0-2)[2i0-3)-c+6l]  =0,  with 
roots  p  =  0>1,2  and  p^=  -lc/2*  Thu8»  aa  in  (2-45) >  we  have: 

♦3(t)  =  (t-i)"^^/2  q(t-l)  (2-95) 


(2-97) 


We  Gsm  also  write  (2*94)  in  the  following  form: 

=  e  !  V  1  '  KfO>i  =  1.2. . 

where i 

Fo(n)=fo(ft*iG/2)=2in(n-ie/2) (n*l*ie/2) (n*2-ic/2)  (2“98) 

P^(n)=f^(n«ic/2«l)=(n-l-ic/2)  (n-2“ic/2)  E  (n»l»ie/2)  (n-  ^(54+313) )  + 

+c(a+i)«.v(v+l)  ]  (2*99) 

P2(n)=f2(n-le/2-2)=-(n-2-ic/2)C(a+b) (n-ic/2-2)^+2b(l-ia) • 

• (n*2-ic/2)*(a+b)v(v+l) ]  (2-106) 

F5(n)sf^(n*lc/2*3)sab(n*3*v*ic/2) (n-2+v-lc/2)  *  (2-101) 

Prom  (2*47) » (2*9$) , (2*96)  we  obtain* 

Wj(y)  =  iIj  (t*!)”^^^^^  q(t*i)  dt  , 

where  q(t-i)  has  a  regular  singularity  at  tsO.  Putting  teiz  we  get: 

V3(y)=  i’^®/^(z-l)"^®/^q[i(z*l)]  dz  ,  Rey>0,  (2*102) 

where  the  mapping  of  the  path  in  figure  (2*2)  i  is  shOiiea  in 
figure  (2*3) •  The  dotted  lines  are  Uie  branch  lines  Of  q[i(z*l)]* 
We  put* 

,-ic/2  «Trc/4  » 

n^jL  n-^o 

where  * 


2“-30 

Then,  according  to  references  1?  (Chapti  XV)  and  18  (Ghapti  III), 
after  substitution  In  (2“102)  and  tenn  by  term  integration,  we 
obtain  for  Rey>  0: 


Wj(y)=e 


( * 


«e 


r(y) 


■[  isinh^ 


Tf  c  e 


.5  r(y)r(n+i«»ic/2) 

IZ  - 

n=o  “  r(y+nfi«>ic/2) 


2  # 

r(n+l^le/2)r(y+l“ic/2) 


1  Z  ‘ 

n^o 


r(y+n+l-lc/2)r(l-ic/2) 


The  expression  in  brackets  is  a  constant  coefficient.  Referring 
to  (2»54) i (2“56)  we  notice  that  as  v^(y)=Wj(y)/r(y)  we  ehoset 


(y)  =  1^  Q  (y)/ny+l* 

g(3) 


n(y)=l* 


i 


,  where: 
(3) 


h=6  “ 


nt 


f{nH*tc/2)f  (y^l^lc/2)  .•ihir(2-^l*ie/2) 

r(y+n+l^ic/2)r(l*lc/2)  r(l»ic/a) (y+l*lc/2) 
(*l)%_r(3*l*lc/2) 


•f. . . . . . 


the  series  being  convergent  for 
that: 


45) ,  (.i)x 


0.  the  equation  also  shows 


(2-105) 


Substituting  (2-105)  in  (2-97)  we  obtain: 


3 

ZT. 


^  "’m'"' ‘^-m  ~  ^  ^  ^  —  0  ,  jiRl, 2,3 > •  •  •  •  ,  (2“106) 


2*31 

(2-107) 


Ti(n)=-1  7”^^  ^ - 7^=  -i(n-l-ic/2)Cn-  f(5a+3b)]+ 

(n-l-lc/2)  (n-2-le/2)  2 


+c(l-ia)^iv(v‘fl) 


P«(n) 


n-2- 


(a*+^)  (n-f2-^ 


tj(n)=lFj(n)=iab(n«3*v«ic/2)  (n^2+v»lc/2)  = 


(&-2-ic/2)  (n‘-3-ic/2)-v(v+l)  ]  ^ 


(2-108) 


(2-109) 


(2-110) 


Next*  we  are  golr^  to  prove  tiiat,  when  a»  b  are  real; 
®n  ^  ®n  ’ 

80  that,  looking  at  (2-56) ,(2-57)  we  should  also  have; 


0  ^ 


(2-111) 


(2-112) 


Furthermore,  since  r(y+l-lc/2)sr(y+l+lc/2)  and 


f  ■W(y^-iyp7f  -iw(y  +lyi)/2  y 

=:e  =  (-l)y,  looking  at  (2- 

and  (2-55)  we  see  that  also  we  have; 


v^(y)  =  Vj(y)  .  (2-113) 

This  relation  holds  not  only  for  Key>  Q  but  also  for  any  y, 
because  In  (2-50),  which  provides  the  analytic  continuation  of 
v^(y)  and  v^(y)  In  the  plane  Rey<Q,  the  coefficients  Pjj(y)  (nsO, 
1»2,3,4)  are  polynomials  in  y  with  real  coefficients. 

In  order  to  prove  (2-111)  we  look  at  the  differential  equa- 


tlon  (2-42)  satisfied  by  around  t^=i  and  t^(t)  around  t^^s 

m-ii  It  IS  an  equation  with  real  polynomial  coefficients  and  its 
two  singular  points  t^“i  and  t^=«i  are  complex  conjugates  with 
exponents  g^=-lc/2=f^.  It  is  then  easy  to  see  that  t^(t)  =  ^ 

which  means  that  C&ee  (2-45)  and  (2*46)]^  Following  the 

Steps  leading  to  (2-165)  for  we  would  find 

.  r(n+l+ic/2) 

•  gQ  tiiat*  by  c^parlson  with  (2-165)  t  relation - 

r(iH^ic/2) 

(2-111)  follows  at  once. 

A  faster  but  not  as  rigorous  way^  is  to  observe  that  (2-112) 
and  (2-113)  hold  at  least  for  ys^^ »  as  equations  (2-54)  and 


(2-57)  show;  since  the  difference  equation  (2-5)  has  real  poly¬ 
nomial  coefficients,  the  relations  should  hold  for  all  y^ 
Turning  to  (2-166)  and  (2-116)  we  can  finds  1  and 


(1)  (-c/2) Cl-i(5a^3b)/2]*c-lac+iv(v^l) 


Ei 


Thus*  the  corresponding  series  expansion  in  (2-88)  can  be  written* 


,-freA+lSglc/2j-i^  ic(a^3b)^i4v(vh|)42C  ^ 


ns2  (1?) 


(2-114) 


Equation  (1-94)  defines  R^(z)  as  follows: 

/X  iz  ic/2r  /C  (c/2i)  (c/21-H)^v(v+l),cb  . 

R  (z)rve^V®/^[l+(§  +  - - )i  .i 

?  e  ^-1  z 


a  e^^  z 


lc(a+3b)+i4v(v<fl)+2e 


i;  .  3  , 


that  is,  the  expression  (2-114)  is  simply:  Ce”  R^(z).  An 
analogous  statement  holds  for  Rj|^(z);  we  finally  conclude  that 
(2-88)  can  be  written: 


i  0<(><fT 


2*33 


R%  (z)  ^ 

2  Izj^oo 


where: 


^4  ^4(?) 


Aj  Rj(z)  *  *Tr<(^<0 

Aj  RjCz)  +  Aj^  R^Cz)  *  ^  =  0  i 


(2*115) 


A^^  =  S  *  ^3  -  ^  *  (2*116) 

For  R2^(z)  we  use  e=v+l  throughout*  for  R2(z)  (when  2v+l  Is  not 
an  integer)  am»v  throughout* 

Since  for  real  a*  b,  z*  R^(z)  are  real  and  R^(z)=  l^(z),  we 

^  2 
mtust  have  A^=  A^^,  or  from  (2*116)  : 

G  =  D  J  for  real  a,  b  ,  (2*117) 


a  fact  that  will  also  be  verified  in  the  process  of  detexmilning 
these  coefficients. 


DETmiNAflON  OF  tHE  COEFTiClENtS  OF  THE  LINEAR  RELATIONS 

In  view  of  (2*116),  our  problem  IS  the  evaluation  of  the 
Coefficients  C  and  P  Of  (2*69).  Referring  to  (2-TO)-(2*75)  we 
have : 

Av^(*3+c)*BV2(-3+c)+Cv^(«3+c)+Pv^(*3+c)  =  0 
Av^(*2+c)*BV2(*2+o)+Cv^(*2+c)+Pv^(»2+cr)  =  0 
AV3^(*l+cr)+iV2(*l+c)+Cv^(»l+cr)+Pv^(*l+o)  =  0 
Avj^(0)+Bv2(o)+Cvj(cr)*pv^(o)  =  1  . 

Solving  for  Af  B,  C,  P  we  obtain: 


(4) 


(s=l,S,3,4)  are  the  cofacters  of  the  last  row  Of  the 


«h«re  njjy, 
detennlhant  P(y)  divided  hy  P(y)  and  where: 


v^(y) 

V 


VjCy) 


(a-119) 


P(y)  is  the  Casoratl’s  determinant  (17  Ghapt*  XII,  18  Ghapt*  I) 
for  the  four  parttoular  solutions  v^(y)  (8=1, 2^3,4)  of  (2*5) •  It 
satisfies  the  following  difference  equation  or  the  first  order 
( Hepiann '  s  the  or  em )  t 


?(y*l)  PQ(y)  _ I _ 

P(y)  ^4(y)  ab(y<*3*v)  (y^4*v) 


(2*120) 


whose  solution  is  given  by  (17  pp.  327*328) * 


P(y) 


(i/ab)y 


^iy) 

r(y+3*v)  r(y+4#v) 


(2*121) 


where  w(y)  is,  in  general,  an  arbitrary  periodic  fimctlon  of  y 
with  period  !•  It  can  be  verified  immediately  that  (2*121) 
satisfies  (2*120).  Since  we  know  how  v_(y)  (8=1, 2, 3, 4)  behave  as 
yp^oo  we  can  find  the  form  of  w(y)  corresponding  to  the  parti¬ 
cular  set  formed  by  these  four  solutions.  It  turns  out  to  be  a 
constant.  As  y^^^  we  use  the  es^ressions  (2*54) *(2*61)  into 
(2-119)  to  obtain: 


(l/ab)^ 


y^r(y+i*lc/2)r(y+l+ic/2) 


.  i/(y+l-ic/2)  “l/(y+l+ic/2) 

(-l/a)^  6  0 

(“i/a)^  (-1A)5  0  0 

In  the  determinant  subtract  tbe  last  column  from  the  third  and 
use  the  difference  as  the  new  third  column#  without  changing  the 
value  of  the  determinant*  The  third  eolumn  becdmea: 

0 

M(y»l) _ 

(y+l«ie/2) (y+l+ie/2) 

0 

0 

so  that: 


(l/ab)^^ 


21  ^ 

y+2^ic/2)  r(y+2+ic/2) 


1  1  0  1 
-l/a  »l/b  1  •»i/(y+i+ic/2) 
(*l/a)_^  (*i/b)^  p  0 
(^1/a)^  (»l/b)^  0  0 


s(i/ab)y*5  £l£  1  r(y»3^v)r(y-»4»v) 

ny+3-v)r(y^4+v)  r(y+2-ic/2)  r(y+2+lc/2) 
But  (7  Chapt.  VIII,  17  pp.  254-255): 


^(y+'f) 

Hy+p) 


,T-p 


(2-122) 


and  this  holds  as  y^oo  in  the  sector  -TT+e<argy  ^n-6.  Applying 
this  to  the  last  expression,  we  obtain: 

1  r(y+3-v)r(y+4+v)  ^ 

P  r(y+2-ic/2)ny+2-^ic/^)  ^  ^  * 


80  that} 


(l/ah)y 


™P  (ah)“'  r(y+3-v)r(y(‘4-Hv) 

Gomparlng  with  (2*lgl)  we  see  that  §(y)  = 
So  finally: 

p{y)  =  (l/ah)y*5  -tig  - -  ^ 

ny+3“v)  rcy+^+v) 


t  a  constants 


We  next  define  the  ao-called  multipliers  (I7  pp^  372*374) 
Ng(y)  (s=lj2i3i4)  of  the  solutions  v  (y)  of  (2*5)  by* 

Ng(y)  =  •^s(y+l)/P4^y^  * 

Here  p^(y)  =  i  and  from  (2*120)  we  obtain*  l/p4(y)=P(y+l)/p(y) ; 
therefore : 

Ng(y)  =  P(y+1)/P(y)  . 

Remembering  the  definition  of  Ug^'y)  (immediately  after  (2*11$)), 
we  see  that  Ug^y^ij  P(y+1)  are  simply  the  cofactors  of  the  last 
row  of  Casorati's  determinant  P(y+l),  We  have  explicitly* 

^  V2(y+1)  V3(y+1)  v^(y+l) 

N^(y)  =  *  ^  V2(y+2)  V3(y+2)  V4(y<^2)  (2*124) 

^  y2(y+3)  V3(y+3)  y^(y+3) 

^  y^Cy+i)  V4(y*i) 

v^(y^2)  V3(y+2)  V^(yf2)  (2-125) 

yi(y+3)  v3(y+3)  y4(y+3) 

^  Vi(y+i)  Vgly+i)  V4(y4'i) 

N3(y)  =  ,  ^  Vj(y*2)  V2(y+2)  V4(y*2)  (2-126) 

'  Vl(y+3)  V2(y+3)  y4(y+3) 


(2»127) 


In  tiie  present  case,  where  p.(y)sl#  N^(y)  (8=:1,2*3»4)  are  simply 

-Q  B 

the  negatives  of  the  cofactors  of  the  first  row  of  P(y)  divided 
hy  P(y)»  The  original  definition  is: 

^sty)  “  P4(y-3.|Ng^y"l)  =  ah(y+2*v)  (y+3+v)Ng(y-l)  i  (2-128)' 

Using  (2-118)  we  can  express  A,  B,  C,  D  as  follows: 


A  s  Cab(y+3-v) (y+4+v)N^(y) (2-129) 

B  =  [ah(y+3-v)  (y+4+v)Ng(y)]y^^2^ifcr  (2-136) 

C  =  Cafe(y+3*v) (y+4+v)N2(y)3^^^^^^  (2-131) 

P  -  Cab(y+3«v)  (y+4+v)N^(y)]y^^^^j^  (2-132) 


For  the  expansion  of  R^(z)  we  use  cr=v+i,  for  that  of  Rg(t)»  when 
Bv-^l  is  not  an  integer,  we  use  o’=-v. 

wording  as  before  for  F(y)  we  can  find  how  Ng(y)  (ss:l,2»3»4) 
behave  as  yj|^^  •  Using  (2-54) -( 2-61) ,  (2-122)  and  (2-123)  we 

find  from  (2-124) : 

r(y+3-v)  r(y+4+v) (ab)^*^ (-l/b)^*^ 

N,  (y)  * 

^  ^rh?"*®  21cyf(y+2-lc/2)r(y+2+lc/2) 

1  1  1 

.  I  -lA*  i/(y+2-lc/2)  -l/(y*2+le/2)  , 

(-lA)^  0  0 


or,  subtracting  the  last  column  In  the  determinant  from  the 
second  cmd  using  it  as  l^e  new  second  column: 
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r(yH>3-v)r(yH-4»v)(ab)y‘*'^(-lA)y*^ 

cr(y+3“ie/2)r(y+3+ie/2) 


,10  1 

»  -l/b  1  -i/(y+2+lc/2)  » 

i/b^  0  0 


or: 


■160 


Similarly  starting  from  (2*125)  we  find  thatt 

Starting  from  (3*136)  we  obtain? 


r(y*3*v)r(y+4‘fv)  (ab)^  ^(«i/ab) 


2lcy^r(y+3^:^ic/3) 


(2*133) 


1 

*l/a 

1/a^ 


(2*134) 


1  1 

0 


and  by  application  of  (2*122)  we  can  find  the  following  two, 
asymptotically  equivalent,  expressions t 


N^(y) 

3  ^  w 


rJ  - 


.v)y’'-2- 

(-1)5' 


r(y'«'4+v)y 


-v-3* 


(2-135) 


Similarly  •' 

N4(y)yi:^  -  ^r(y+3*v)y^”^*^^/^A;-  ^r(y+4+v)y*^‘‘^’*‘^°/2.  (2»136) 


The_ Ad.lolnt  Difference  Equation:  It  Is  a  well-lmown  fact 
that  the  multipliers  N_(y)  (s=l,3,3,4)  are  Independent  solutions 

3 

of  the  difference  equation  adjoint  to  (2”5),  i*e*  they  satisfy 
(17  pp.  372*374,  18  Chapt.  I) * 


2«39 

4 

%  Pjjj(y+4^m)N(y+4«m)  =  0  * 

ffls© 

dr: 

<l4(y)N(y‘^4)+q3(y)N(y+3)*^q2(y)N(y^2)+q^(y)N(y4l)  + 


+qo(y)N(y)  =  0  ,  (2*137) 

where : 

Q4(y)=PQ(y+4)  =  1  (2*138) 

q3(y)=Pi(y+3)  =  2a  (2*139) 

q2fy)-p2(y+2)^(y‘^4)  (y+3)+a^“v(v+l)  = 

=  (y‘»-2)(y+3)+2(y+2)+a^+2*v(v+l)  (2-140) 

qi(y)=q?3(y+l)=(a+b)  (y+l)  (y+2)+(5a+3b)  (y+l)+3c^ 

+(a+h)£6*v(v-t*l)i  (2-141) 


9o(y)sp^(y)*ah(y+3*v)  (y+4+v)=ab[y(y+l)if6y*l2*v(v+l)i  .  (2*142) 

Applying  the  methdd  of  Laplace's  transformation  as  before »  we 
form  the  functions: 

^gd)  =  t^+(a+h)t+ah  s  (t+a)(t+b) 

<(•1(1)  =  2t^+(5a+3b)t*6ab 

(^g(t)  =  t^+2at^+i:a^+2*v(y+l)3t^^C3c+(a+b)  (6*v(v+l) )  3t+ab[l2*v(v+l)  ] 

The  corresponding  differential  equation  for  ilr(t)  Is: 

(t*a)(t-i»b)t%"(t)*4j(t)tt'(t)+^Q(t)^f(t)  =0  (2*143) 

and  has  regular  singula*  points  at  t=0,  t=*a,  t=»b  wid  an  Irregu*^ 
iar  singularity  at  t?=®^*  Lndlelal  equation  around  t=0: 


2«4d 


tl3p(p«l)»6afep+:ab[l2«v(v+l) ]  =  6  with  roots  p|=  p|=3-v. 

Slhce  V  is  positive i  for  Rey>  v»3  the  GOrrespOhdihg 

is  zero  for  any  solution  of  (2-143) •  the  indiclal  equation  around 
is: 

,  ,  ^  ^  2h^«(5a+3fe)h+6ab 

(-#)/(lDc)  ]g=d  with  roots  0  and  P«=l-  :  — v-:-~-r-.r--Tr7  s 

2  no 

=  0.  There  exists  one  solution  tg(ti=  l<fd|(t-*'p>>fdg(t+p)^+* . . . . 
of  (2-143)  andi  correspondingly,  a  solution  of  (2-13?)  in  ■tile  form: 


analytic  for 


v-3  . 


The  path  £  is  shovm  in  figure  (2-4)^  As  in  (2-20) ^  it  can  easily 
be  shown  that  l('^2»‘t)  0.  According  to  (2-52)  *  (2-53) » 

NbrlundU  theorem  and  the  statement  on  page  2-20 >  we  can  expand 
the  above  expression,  as  *  as  follows: 


r(y) 

r(y%l) 


Comparing  with  (2-134)  we  see  that  the  above  solutipn  is  propor¬ 
tional  to  ^^(y). 

Similarly,  the  other  solution,  corresponding  to  the  singular 
point  t=-a  of  (2-143) »  would  be  identified  as  proportional  to 
Nu^(y),  We  are  Interested,  however,  only  in  the  solutions  N^(y) 
and  N^(y)  of  (2«?i37) ,  on  which  0  and  D  depend.  The  coefficients 
At  B  do  not  appear  in  the  asymptotic  expansions  of  R^(z)t  R2(z), 
We  will  find  N^(y)  and  ^^^(y)  explicitly*  We  putt 


N(y)  ^  r(y+4*v)  M(y)  .  (2-144) 

Substituting  in  (2-137)  we  obtain: 


Q4(y)M(y+4)+Qj(y)M(y+3)'*^(7)M(y+2)*t%(y)M(y4l)+ 


+Qo(y)M(y)  =  0  ,  (2»145) 

where,  with  the  use  of  ( 2‘-138)*( 2*142)  *  we  have: 

Q4(y)-q4(y)  ny+8+v)/f(y+5+v)  =  (y+4)  (y+s)  (y+6)‘:^3(v+l)  (y+4)  (y+5)  + 
+3(v+l)  (v+2)  (y+4)-f  (v+1)  (v+2)  (v+3)  (2*146) 


Q2  (y )  -^2  ( y )  r(  y + 6+ v )  /  r(yH*5+ V ) = ( y+2 )  ( y +3 )  ( y  ^^4 ) + ( v^3 )  ( y+2 )  ( y+5 )  + 

+ia^+6*v(v«l)l(y+2)+(  V‘f3 )  [  a^+ 2*  v  ( vf  1 )  ]  ( 2*148) 

(y+i)  (y»*2)+(5a+3b)  (y+l)‘t-3c+ 

+(a+b)[6*v(v+i)]  (2»i49) 

QQ(y)saby+ab(3*v)  .  (2*150) 

Applying  the  method  of  Laplace's  transformation  we  form  the 
functions : 

4^=(t)=t\t^=t^(t2+l)  (2-151) 

^2 ( t ) =3 ( v+1 ) t^+2at^+ ( v+3 ) t^+ ( a+b ) t  ( 2*152) 

^^(t)=3(v+i)  (v+2)t^+4a(v+2)  t^4-[a^+6«»v(v*l)  it^+ 

+(3a+3b)t+ab  (2-153) 

^p(t)=(vti) (y+a) (y+3) t^+2a(v+2) (y+3) t^+(y+3) [a^+2*v(v+l) ]t"+ 

♦C3c+(a+b) (6-v(v+l) ) ]t+ab(3*y)  (2-154) 

The  corresponding  differential  equation  for  ^(t)  and  expression 
I(y»t)  are: 


t2(Al)tV*(1s)-4g(t)tV(t)+4j(tH<l''(t).fj,(t)*(t)=6  ,  (4-15S) 

m=d  dt  m=d  at 


+\jr'^(t)  =  tCt)  (!)<♦»  (y+i)(j)g(t)+t<fg(t)  + 

(y+2)  (y+l)  4^  ( t)  +2(y+2)  tf^( t)  ( t)  ] t)  + 

+  (y+2)^3(t)+t<|i|(t)]+t''(t)ty‘^^f^(t)  .  (2-156) 


Eqiiatldn  (2*155)  Ms  three  regular  Bingularities  at  t^=ii  t^=-l, 
t=6o  and  an  irregular  singularity  at  t=0»  For  Its  solutions 


»e“'  ,  0^=1  and  obtains 

n=o 


around  t=d  we  puts 

♦  '(t)se“‘/‘|  (2.1S7) 

nso  ^  n=o  “ 

♦"(t)=e®/^  JE  (n^3)  (n^t^^-Do  2  (n+p-l)ct^t®'*’^“^+ 

n=o  nso  ® 

+m^e®/^  1  ( 2-158) 


n=o 


V"(t)=e®/^C  2  (n+P)  (n+P-1)  (n+p»2)a„t^*^"^-3m  1  (n+p-1)  (n+p-2)  • 
n=o  “  nso 

•  a  2  (n+P-a) o  2  a  .  (2-159) 

"  n=o  -  h=o  ” 

Substituting  into  (2-155)  and  eliminating  the  common  factor  e®^^t 
we  obtains 

(t^+l)[  2  (n+p)  (n+P-1)  (n-*-P-2)o  t®’'''«3m  2  (n+P-1)  (n+P-2)a 

n=o  ”  n=o  “ 

*3rar^  (n+P-2)o  t®-m^  2  t^+2at~+(v+3)  t+(a+b)  3 • 

n=o  n=0 
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•I  ^  (a+p)  2  (a+#*i)a^t*^+m^  ^  o  t^'^^3+f5(v+l)  ‘ 

a=o  nso  “  a=6  “  ^ 

•  (v+2)  t^+4a(vs^2)  t^+Ca^+6«v(v«l)  3t^+(Ba‘f3te)  t+aTs^  •  [  %  (n+P)Oj.t^« 

'J  a=o 

*m^  (v+2)  (v+3)  t^+2a(v+3)  (v+3)  t^+{v+3)  [a^+2“v(v+l)  ]• 

a=6  ^ 

•  ACBe+Ca+b)  (6-v(v+l))]t+ab(3«v)?  2  =  0  « 

*'  ri=  o  “ 

The  lowest  power  of  t  la  this  expreBBloa  Is  t~^*  Equatlag  the 
coeffleieatB  of  t“^  aad  t®  to  O  we  obtaia! 

t"^:  «»m^a^-(a+b)m^a^‘-abma^@  0  with  roots  a|=Oj  mg=«b  * 

t®  t  «m^aj^+3m^(0^2)aQ-(a+b)m^a^*(v>f3)ffl^o^+2in(a+b)  ($-l)oQ»abmaj^- 
•*»(5a+3b)maQ+ab$a^*ab(3‘S’v)a^  =  0  ,  or 
(3$‘^v-9)m2+[2(a+b)(:$^l)-5a-3b3ffl+abO+v-3)  »  0  * 

For  Oj^sO  we  get:  3»v  . 

For  m2=«a  :  {J)^^^v^9)B,*2{3.*'b)  or  Pg®  v+a. 

For  m^=-»h  s  (3P2’^v*9)b«'2(a+b)  (p^‘^l)*5a^^3b+a(Pj+v^3)sO»  or  P^s  v+4. 
Thus,  three  normal  asymptotic  solutions  around  t^  are  obtained: 

a=l 

nsi 

♦iii't’t^o  t’'**  Ci+z  <■'“*’ t”]  . 

n=l 

Upon  any  ray  drawn  from  t^O,  for  which  Re(»a/t)XOt 
Re(»b/t)^0  and  for  suiy  y  we  have: 

llm^^p#i|rg(t)=llm^^Qt5^*H^(t)=:llm^^Qty*%''(t)  =  0,  for  s-|I,III. 


1  ! 

‘  i 
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If  we  further  restrlet  y  to  values  for  which  Re(y+3-v)>  0,  or 
ley  >  v-3 ,  we  are  going  to  have : 

llin^^Qtyt|(t)=llm^^Qty*^Hj(t)=lia^_^Qty*2^j(t)  =  0  . 

fherefore,  for  any  path  i  in  the  t-plane,  starting  and  ending  at 
t=o  along  rays  satisfying  the  conditions  Re(“a/t)<ro,  Re(-b/t)^Oi 
while  Rey>v«3i  we  are  golr^  to  have  l(t,t)|^  s  0  for  any  solution 
^(t)  of  (2-155).  This  is  easily  verified  by  looking  at  (2*156)* 
(2*151) -(2-1S4)  .  furthermore,  the  sajie  resist  is  still  valid,  if  _ 
the  path  approaGhes  t=o  along  either  of  the  two  limiting  rays: 
Re(-a/t)=0*  Re(*b/t)4  0  or  Re(*a/t)^0,  Re(«b/t)=0*  if  y  is 
restrleted  by  Rey > max(v*3,-v*2) *  a  condition  that  reduces  to  the 
previous  one,  Rey>v*3*  if  v^l/2.  Again*  a  look  at  (2-151)* 
(2-154),  (2-156)  and  (2-15?) *(2*159)  verifies  this  statOTent 
immediately.  As  before,  these  observations  will  serve  to  fix  the 
path  of  Integration  in  the  t*plane. 

The  indlclal  equation  of  (2*155)  around  tj|^=i  is: 

t^<|)o(t) 

P(g-l)  ( $-2) *[ ]  p ( p-l)  =  0  with  roots  0*  1  and 
t^(t#i) 

^4- 

indlclal  equation  around  t,=*l:  5 ( p*l)  ( 5*2) .  ,$(^-l)  = 

^  t^(t-l) 

~  0  with  roots  0,  1  and  0^^  v+2+ie/2. 

There  exist  two  solutions  of  (2*155)  in  the  form: 
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'(3) 


(y)  -  e 


r(y+4+v) 

’  2frl 


at 


(2-163) 


whe?e  the  posslhle  paths  of  Integratioh  are  shown  in  figure  (2-S) 
fhe  so  defined  funotions  are  analytic  for  Rey>v«3  and  as 

in  figure  (2-5)  I  or  for  Rey>  max(v«3i-v«2)  and  and/or 

tangent  to  the  liffiiting  rays  as  Wo*  For  v>l/2  both  oases  are 
Equivalent.  As  y|^6d  we  can  expaiia  the  atinve  fimetlona  ao  hofftra 


in  the  following  forms: 

r(y+4+v)r(y) 

^  r(y4Vi»i3»ie/2) 


K. 


(4) 


y+v+3- 


d^i^  r(y+4‘fv)y 


/  4  ^  .  r(y+4+v)r(y) 

“(3)  °<>ny.vt5*l«/8) 


(- 


K 


«Vi.' 

(3) 


y+v+3+ic/2 


+ . 3  rs^ 


A»c^(-i)^  r(y+4^v)y"^"^” 


comparing  with  (2-135) • (2-136)  we  see  that  the  particular 
functions  N^(y)  and  N^(y)  we  are  looking  for*  are  obtained  if  we 
take;  d.s  c  s  ^1/2  ,  i.e. 

r(yM'fv)  1  r  ^  T 

H^(y)  - - - -  pj  (2*16*) 

(y+4+v)  T  r  , 

N^(y)  =  ,  ^  ™  2^  •  (2-165) 

They  are  analytic  functions  of  y  for  Rey>v»3  and  as  in 

figiire  (2-5) »  or  for  Rey >  max(v»3.-v-2)  in  the  limiting  case, 
mentioned  above.  Their  continuation  to  the  left  is  provided  by 
the  difference  equation  (2-137)  itself,  i.e.  by? 

N(y)=-  «  1  L.U  YCN(y*4) »2aN (y*3) ♦q^Cy) W (y*2) » 

ab(y+3-v) (y+4+v)  - 
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+%(y)N(y+l)3  .  (2*166) 

Therefore,  they  have  simple  poles  at  the  points: 

y  =  v«3»n  and  y  =  -v«4*n  *  n  =  0»l*2,**i  ^  (2*167) 


Another  remark  can  he  made  at  this  point!  When  lef|=v+2* 
“Imagc/2>*1  and/or  lep^=v+2H»imagc/2>  *1*  the  integration  around 
the  Small  circles  aurroundlng  the  points  *1  and/or  1,  In  and/or 
yields  0,  as  the  radius  of  the  circles  Is  allowed  to  go  to  0* 
We  can  then  replace  the  corresponding  paths  by  Integrations  along 
the  lines  and/or  as  shown  in  figure  (2*6).  The  only  change 
in  the  expressions  (2*164) , (2*165)  for  N^(y) ,  Nj(y)  amounts  to  an 
appropriate  constant  factor  In  front,  easily  determined  in  each 
case. 

Explicit  solutions  of  (2*137)  for  Nj(y)  and  N^(y)  can  he 
found  hy  the  method  used  previously  to  obtain  explicitly  the 
solution  v^(y)  of  the  difference  ecjuatlon  (2*5).  The  steps  are 
exactly  the  same.  The  factorial  series  obtained  for  Vj(y)  is 
convergent  for  Rey^  o.  in  the  present  case  the  series  obtained 
are  at  least  asymptotic  and  the  process  is  at  least  formal.  We 
start  from  the  defining  Integral  representations  (2*164) , (2*165) 
and  proceed  to  obtain  explicitly  the  series  developments  (2*160) 
and  (2*161)  for  ^^(t*i)  and  t^(t+i),  which  satisfy  the  differen* 
tial  equation  (2*155)  with  coefficients  ^g(t)  given  by  (2*151)* 
(2*154),  The  final  results  are  as  follows: 


N4(y) 

3 


r(y+4+v)r(y) 

- li— T(y) 

r(y+v+3+ie/2) 


T(y)  =  1+ 


d. 


y+v+3+ic/2  (y+v+35ic/2) (y+4+v#lc/2) 


t.. 


(2*168) 

(2*169) 


Reemrenee  formula  for  the  coefficients  d^: 


^  ®  .  (2-170) 

X  =  ri#v?iG/2  (2-171) 

f]^(x)=llx^'*[14v«*‘17?l(7a-3fe)  ]x+4v^+8v-ac?l(4av+3d)  (2-172) 

-rg{x)=»2Sx^4‘C39v+102fl<3^l7^)  Ix^^C  lSv^*81v*89*^«3a^|i 

+39) )  ]x+2v^+l0v^+(l2«a^-at))v«3ac?i[  (Sa+fe)  v^+(lla‘»»b)  v+3e]  (2-173) 

T'3(x)=xi30x^-{56v+228+l(b-19a)]x^+{3(llv^‘f85v+l72-a^)|ii(3b- 

-a(12v+50) )  ix-6v^-65v^+(2a^-24S)  v+9a^-330fiE  (7a+b)  v^+(55a-fb)  v+ 


+I08a-8b]^  (2-174) 
( x) =x ( x-1 ) ( X- V -5 ) C -20x^+ ( 24v+l22?ll0a) x+a^-Tv^-STv- 

-174+12a(3v+16) ]  (2-17S) 
‘*'5(x)=x(x»l)  (x«2)  (x-v-6)  (x-v-5)  (7x-4v-25j:21a)  (2-176) 
t'g(x)s-x(x-l)  (x-2)  (x*3)  (x-v-7)  (x-v-6)  (x^v-5)  .  (2-177) 


Throughout  the  above  relatlonB  the  upper  sign  is  used  for  N^(y) , 
the  lower  for  N^(y). 

Differentf  but  equivalent,  expansions  for  N^(y) ,  N^(y)  can 
be  found  by  employing  Boole’s  operational  method  for  solving 
difference  equations  with  rational  coefficients,  as  explained  in 
reference  17  (Chapt*  XIV  pp.  434-461),  it  Is  essentially  the  me¬ 
thod  of  undetermined  coefficients,  completely  analogous  to  the 
method  of  Probenius  for  differential  equations.  The  method  is 
tedius  and  lengthy  but  straightforward.  The  final  results  are  as 
follows : 
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N4(y)  =  “ 

3 


(|i)^  r(y‘*‘4‘^v)r(y+3*v) 
2  r(y+65le/2) 


=  1+ 


y+6+iG/2  (y+6lie/2)  (y^i-T+ic/g) 


+  » 


( 2*178) 

(2*179) 


Recurjfenee  formula  for  the  coefficients  d  i 

n 


^o"^  ^  *»  0  »  J=lr2i3,4*44  (2*180) 

X  =  2»v*n+ic/2  (2*181) 

fg(x)=x(x-l)  (x*2)  (x+2v+l)  (m2v)  (x^2v*l)  (2*182) 

f|(x)=2(x*l) (x*2) (x+2v) (x+2v*l) (3x+3v-3?la)  (2-183) 

f 2 (x) = ( x-2) ( X+2V-1 ) C 14x2+4 ( 7v*10ll2a) x+12v^*42v+30-a^l 

+14a(2v*3)]  (2*184) 


f^(x)=16x^+C48v*90*l(lla*b) ]x^+2[20v^-94v+8|*a^li( (ila*b) (v-2)+ 
♦b) 3x+8v^*84v^+(184-2a^) v+4a^-l08?lC  8av^+(7b*47a) v+ 


+44a*8b] 

f ^(x) =9x^+[ 18v-39+i2 (3a-blx+8v^^40v+42*ac+l[ 2( 3a*b) v+ 

+7b-15a3 


(2-185) 

(2-186) 


The  si^  convention  is  the  same* 

The  inverse  factorial  series  In  all  these  expressions  are  at 
least  asymptotic  for  large  |y|  *  It  is  also  obvious  that  with  real 

a,  b;  Then:  N^(y)  =  N^(y),  or  for  real  y:  N5(y)=f^(y) 

Referring  to  (2*131) » (2-132)  and  with  y=a»4=v«3  for  Rj_(z),  y=-v-4 
for  R2(z)  (when  2v+l  is  not  an  integer),  both  of  which  are  real 
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values,  we  Gonclude  tKiat!  C  =  5  ^  a  fact  which  was  Inferred  pre¬ 
viously,  equation  (2-117)*  through  different  considerations* 

It  may  also  be  pointed  out  that  different  expansions  for 
N/,.(y)  *  Nj(y)  ean  be  obtained,  which  may  prove  better,  from  the 
computational  point  of  view,  depending  on  the  particular  values 
Of  a  and  b  imder  consideration,  if  in  (2-137)  we  use  the  general 
substitution: 


N(y)  =  rCz+h)  M(z)  ,  z  =  y+r  (2-187) 

and  solve  for  M(z),  employing  either  of  the  above  methods*  The 
parameters  h  and  r  can  be  chosen  conveniently  in  each  case  so  as 
to  optimize,  from  the  computational  point  of  view,  the  expansions 
obtained*  loth  expansions,  given  explicitly  above,  were  used  for 
the  computations  in  chapter  3,  Pi^T  I,  yielding  the  same  values 
for  G  and  p.  The  fomer  proved  better  in  certain  cases*  others, 
based  on  the  change  of  variable  (2-187),  were  also  used  for 
comparison  and  cheek  on  the  results* 

The  evaluation  of  c  and  p  is  based  on  equations  (2-131), 
(2-132),  with  y=o’-4sv-3  for  R^(z),  y=or-4s-v-4  for  R2(z),  when, 
in  the  latter  case,  2v+l  is  not  equal  to  a  positive  integer*  In 
general,  these  values  will  not  be  useful  for  dlreet  evaluation  of 
N(v-3)  and  N(-v-4).  However,  we  can  evaluate  N^(y) ,  Nj(y)  at  y, 
y^+l,  y^2,  y^3,  where  y  is  adequately  large  and  then  use  the 
difference  equation  (2-166)  itself  to  obtain  values  for  N(y-l), 
N(y-2)  etc*,  up  to  n(v-3)  and  N(-v-4)* 

All  these  considerations  c^  be  expressed  in  another  more 
compact  and  general  form,  which  will  also  prove  necessary  in  the 
next  section,  where  the  case  of  integral  values  for  2v+i  is 
investigated* 

Instead  of  the  initial  conditions  ( 2-70) -( 2-73 ) ,  we  make  use 
of  the  general  ones  given  in  (2-68)  •  v(n+o)  s:  a^^  for  all  integers 
n*  It  has  been  proved  that  they  are  equivalent  to  the  four 
expressed  by  (2-70) -(2-73) •  The  function  v(y)  is  defined  in  (2-69) 


2«S0 

Starting  from  (2«119),  tne  definition  of  Casoratl's  determinant^ 
we  multiply  the  first,  second,  third  and  fourth  columns  hy  A,  B, 
G,  D,  respectively,  add  and  use  the  sum  as  the  new  fourth  column. 
Applying  (2»69)  we  obtain: 


This,  Incidentally,  shows  that  the  constant  D  is  equal  to 
the  constant  ratio  of  two  of  Casoratl's  determiinants  of  the 
difference  equation  (2-5) :  D  =  »  where  P(y)  corresponds 

to  the  four  particular  and  independent  solutions  Vg(y),  s=l,2,3,4 
of  (2^5),  while  P^(y)  corresponds  to  the  set  of  solutions  v  (y) , 
s=l,2,5  and  v(y)  of  (2*5).  From  (2*188)  we  obtaint 


Ds 


^2(7)  v^(y) 

V2(y4l)  V^(y+1) 

V2(y+2)  v^(y+2) 


7|(y)  V2(y)  v^(y) 

v(y+l) 

Vi(y)  ^2^7)  Vj(y) 

V3^(y+1)  V2(y+1)  v^(y^l) 

Vj^(y+2)  Vg(y+2)  Vj(y*2) 

V3_(y+5)  V2(y+3)  v^(y+3) 

P(y) 

Vj^(y*3)  V2(y^3)  Vj(y»-3) 

v(y) 

P(y) 


v^(y+i)  V2(y+1)  V^(y*l) 
v^(y+2)  V2(y+2)  v^(y+2) 
Vj(y+3)  V2(y+3)  v^(y*-3) 


(2*189) 


We  call  (a),  (g),  (y) ,  (b)  the  four  terms  of  (2*189).  Making  use 
of  (2*177)  and  (2*120),  i.e.  P(y+l)  =  P(y)/p^(y) ,  we  have* 

(e)s|;y(y+3)/P(y)  ]P(y*l)N^(y*l)=v(y+3)p^(y-l)N^(y-l)  (2*190) 

(6)  =  *V(y)N^(y)  .  (2*191) 

For  (#)  we  observe  that  for  0=1,2, 3  we  have  from  (2-5): 


Vg  (y+3 )  =  •  c  P3  ( y  )  Vg  ( y+^)  +P2  ( y*l)  Vg  ( y+l )  +2aVg  (y )  +Vg  ( y«l )  ]/p2^(  y  *1 ) . 


Then  we  sufeitltute  into  the  last  now  of  the  detefUlnant  in  ($) 
and  hreak  it  up  into  four  determinants  with  the  same  upper  two 
rows^  as  in  (|)>  and  with  last  rows  containing  each  of  the  four 
terns  in  the  ahove  expansion,  respectively.  The  second  and  third 
of  these  determinants  vanish  because  they  have  two  proportional 
rows.  There  remains i  —  - 


O)^ 


v(y+2)  p^ly*!) 

p(y)  p4(y^i) 


vi(y) 

v^(y«fl) 

v^(ye2) 


VgCy)  v^(y) 
VgCy+l)  v^Cy-fi) 
VgCy+a)  v^(y+2) 


Vj^Cy)  VgCy)  v^(y) 

•  V|(y-fi)  VgCy+l)  v^(y+l) 

Vj(y«l)  VgCy^l)  Vj(y-l) 


Interchaniing  first  and  third  rows,  then  second  and  third,  in 
the  last  determinant,  we  finally  obtain? 


(») 


^(7*2)  p^Cy"!) 

P(y)’  P4{y-i) 


f(y*i)N4(y«i)+ 


^(y+3)  1 _ 

P(y)  P4^y*i) 


P(y-2)N4(y-2)= 


“  v(y+2)p3(y*l)N^(y«i)+v(y*2)P4(y-2)N^(y--2)  .  (2»192) 

Similarly  for  (r)  we  substitute  the  elements  of  the  first 
row  using? 


Vg(y)  =  «CP4(y)Vg(y+4)4.p^(y)Vg(y«t.3)+P2(y)Va(y+2)+2aVg(y+l)], 


with  8=1, 2, 3.  We  finally  obtain? 


(y) 


P4(y) 


Vj(y)  Y^iy)  y^(y) 
Vj^(y+2)  V2(y+2)  Vj(y42) 
▼j^(y+3)  V2(y+3)  Y^{7’¥’i) 


v(y+l) 

P(y) 


2a 


Vg(y+1)  v^(y+l) 
VgCy+a)  V|(y+a) 
Vg(y+3)  v^(y+3) 


P4(y)P(y+i)N^(y+i)* 


m  f  ^  4  \ 

^  L£™£aP(y)N^(y)=-f(y+i)N^(y+i)-2ay(y+i)N4(y)  .  (2-193) 

P(y) 

Substituting  (2-190) «•( 2*193)  into  (2*189),  we  finally  obtain: 


^-s-vCy+i)  P4(y-2)N^(y-2)  +[v(y+3)  P4(y*l)  +vCfSlTp^Ty^OTK^(y“i)  - 
“Cv(y)+2ay(y+l)  ]N2^(y) *v(y+l)N4(y+l)  *  (2*194) 

For  y^+o  and  relation  (2*68) :  v(n+(?)=aj^  for  all  n,  we  obtain: 


P  =  aj^^2P4(n*2+o)N2^(n*2%d)+[aj^^3P4(n*l+or)+aj^^2p3(n-l+dr)  ]N4(n*l+o') 

*Ca^+2a  a^^3^]N^(n*o)*aj^^3^N4(n+l4o)  .  (2»19S) 

For  an  appropriately  large  value  of  n,  the  four  values  of 
N4(^)  appearing  in  (2-19S)  can  be  evaluated.  As  a  check,  we 
evaluate  D  using  (2*193)  for  10  or  12  values  of  n,  either  conse¬ 
cutive  pr  not.  For  n**3  we  simply  get:  P=p4(C*4)N^(d*4) ,  as  in 
(2*132).  C  is  given  by  the  same  equation  (2*195) »  if  NjCy)  Is 
substituted  in  place  of  N^Cy). 


2v-:>l  IS  EQUAL  TO  AN  INTEOER 

Asymptotic  Expansion  of  R^(z)  for  Large  Isl  :  For  the  coef* 
ficlents  Bjjj  of  the  logarithmic  solution  ^2(2^)  we  obtained  in 
Ghapter  I,  PART  II,  the  recurrence  formula: 

B/^(m+v+l)+ig^^^f^(m+v)+B^^3f2{m+v,l)+B^^3f3(m+v-2)^„j,4f4(®+v-3) 

=  ajuFQ(m+v+l)+«^_^Fj^(m+v)+ajp^j^2(m«fv-l)  .  (1*55) 


2*53 

The  h^’s  are  coimeeted  to  the  Bg^'s  throiaghl 

n  «  av+l+o  .  fiavfl+B*  ®a  * 

The  initial  conditions  for  these  fO£inulas  are  foxmd  from  (1*53) 
and  (1*54)  to  fee: 

B.gv-A®  *.av.3=  ®-av.g=  0  •  B,gv-1=  ®6  ^ 

With  1^=  fe«.,^T=  0  and  fe^  as  in  (1*51) »  these  conditions  can  fee 

0  2V+iL  O 

replaced  fey  the  equivalent  set: 

®o”  ®  »  ®*1“  ^o^iv  ♦  ®*2^  ^o*^2v*l  * 

^2v«2^  ^o^2v*2  ’  (2*197) 

where  fe^  and  dg^*  ^2v*l*  ^2v*2  definite  nuafeers  defined 

previously*  equation  (1-50)*  It  is  easily  checked  that  (2*197) 
lead  to  the  values  (2*196)  through  the  difference  equation  (1*55) 
and  the  special  value  of  fe^  in  (1*51)* 

We  next  write  m+4  for  m  in  (1*55) •  Reaemfeerlng  from  (2*68) 
and  (2*69)  that  in  this  case  a^gS  v(m+v+l)  we  ofetaln: 

fo(m*v+5)Bgj^^*fj^(m+v+4)ig^^3H.f2(m+v+3)B^^2*f3(m+v*2)Bgj^^+ 

♦f ^( m+v+l)Bjg“  Fq  (m+v+5 ) V  (m+v+5 )  +Fj^  (m<fv+4)  v  (m+v+4)  ♦ 

+F2(m*v+3)v(m+v-f3)  .  (2*198) 

We  introduce  at  this  point  in  place  of  m+v+l  the  general  varlafele 
y  and  in  place  of  the  function  v(y)  such  that* 

v(m+v+l)  =  Bjg  for  all  integral  values  of  m.  (2*199) 


Equation  (2*198)  transforms  Into* 
4  2 
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This  equation  for  y=  m+v+l*  m  being  an  integer,  readily  reduceB 
to  (2*198)*  if  (2*199)  is  aiao  used.  We  write  it  finally  as 
follows : 

P4(y)v(y+4)+pj(y)v(y+3)‘«-Pg(y)v(y+2)H‘Pj^(y)v(y+l)+PQ(y)v(y)  = 

=  qg(y)v(y+4)4‘qj^(y)v(y+3)‘«‘qQ(y)v(y+2)  »  (2*200) 

wbere,  with  the  use  of  (1*38) “(1*43) ,  we  havej  _ 


P4(y)=^o(y‘‘‘4)=ab[  (y+4)  (y+3)*v(v+l)  ]=ab(y+3*v)  (y44#v)  (2*201) 


p3(y)=f^(y+3)  =  (a+b)  (y+3)  (y+2)+c(y«i‘3)*(a<fb)v(v+l)  (2*202) 

P2(y)sf^(y+2)»(y+2)  (y+l)-t-a^*v(v+l)  (2*203) 


P|(y)=f3(y+l)^2a 


(2*204) 


Po(y)sf4(y)s  1 

^pi.(y)’ 

qg(y)^o(y*4)=  ^abC2(y+4)*i]=  *  *|i~ 


(2*205) 

(2-206) 


Ql  (y )  (y+3 )  =-  ( a+b)  [  2  (y+3)  -1]  *e=*2  ( a^b )  ( y+3 )  +26= 


qQ(y)^Fg(y+2)=*2(y+2)+is*2y*3- 


dP3(y) 

Ty 


4pg(y) 

4y 


(2*207) 
( 2*208) 


The  last  expressions  for  qg(y)  as  *[dPg^2(y)/dy]  ,  s^0,l,2  are 


found  from  (2*201) *(2*203)  by  differentiation,  in  (2*200)  the 
fiuiction  v(y)  is  the  solution  of  the  difference  equation  (2*5) 
which  corresponds  to  R^(z),  i,e«  the  one  that  can  be  expressed  as: 


y(y)  Av^(y)+BVg(y)*Cv^(y)+Dv^(y)  (2*209) 

and  which  reduces  for  y;ai+v+i  to  y(n’*-v+i)=a^,  the  coefficients 
of  R2(z):  Vg(y)»  s=1,2>3»4»  in  (2*209)  are  the  preyiousiy  found 
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fouf  IMepeiadent  particular  solutiGns  of  (2*5)‘ 

Gomparlson  of  Pg(y).  s“0»lj2i5»4,  as  given  by  (2«201)« 

(2*205) »  with  (2*6) *(2*10)  Shows  that  they  are  identical  with  the 
Po(y)  of  (2*5) •  fhat  is,  the  corresponding  to  this  case  difference 
equation  (2*200)  is  an  inhomogeneous  difference  equation  whose 
homogeneous  part  is  identical  with  (2*5) *  We  are  looking  for  the 
special  solution  v^(y)  of  this  equation^  whieh  satisfies  the 
.Gfinditlons  (2^499)  ,  l*e.»  -  - — 

Vj^(n+v»t>l)  s  B^=  for  all  Integers  n  i  (2*210) 

As  before,  only  four  of  these  eondltlons  are  sufficient,  i.e. 
from  (2*197)! 


v^(v+l)=0  ,  Vi(v)-B^l=b^d2^  ,  v^(v»l)=B_g=b^d2^^l  , 


More  generally  we  can  use  (2*210)  for  n»  n+1,  n+2,  n+5»  where  n 
is  any  integer  and  ^n+2v^l*  coefficients  appearing  in  the 

definition  (1*57)  for  ^^(z). 

The  general  solution  of  the  Inhomogeneous  difference  equation 
(2*200)  consists  of  the  general  solution  Of  the  homogeneous  equa- 
tioni  i.e.  of  (2*5) f  plus  a  particular  solution  Vg(y)  of  the  in¬ 
homogeneous  equation  (17  Chapt*  XII).  That  is,  we  are  going  to 
have ! 


Vj^(y)  =  iVj(y)+Fv2(y)'fOv^(y)+Hv^(y)+Vj(y)  =  V(y)+yg(y)  ,  (2-211) 


where  E,  F,  G,  H  are  constants,  which  are  going  to  be  determined 
so  that  the  initial  conditions  (2-210)  (Just  fovu*  of  them)  are 
satisfied. 

In  order  to  obtain  a  particular  solution  yg(y)  of  (2-200)  we 
consider  the  function  vCy)  as  defined  in  (2-209) •  It  satisfies 
equation  (2-5) »  i*e. 


P4(y)^(y+^)+P3(y)^(y‘^3)+P2(y)v{y+2)+p^(y)^(y^‘l)+p^(y)^(y)  =  o 


Differentiating  with  respeet  to  y  and  remejaiberlng  from  (2“ 

h 

dy 


,  .  dPiiCy)  ap7;(y) 

(2»208)  that  =«qo(y)  *  -  =- 


dy  '^2 
dPi(y)  tPoly) 

while  -  =  6>  we  obtain: 


,  .  d^(y) 

%(y)-  -  =“tgW. 


dy 
dT(y+4) 


dy 


IT" 


dv{y+5) 


♦PoCy) 


dV(y+2) 

dy 


+Pi(y) 


dv(y+ir 

dy 


■’■po^y)'dy  ‘  (2*212) 

After  eofflparlBon  with  (2-200)  this  equation  shows  that  a  parti* 
cular  solution  of  (2*200)  is: 


(2-215) 


dv(y)  dv,(y)  dv^(y)  dv^(y)  dv^(y) 
Vg(y)=  s  A“^^***«  +B—^  +0***^ -  +D 


dy 


dy 


dy 


where  a,  S,  C,  D  are  the  definite  eonstants  corresponding  to  the 
first  solution  R^(z)  with  osv+l*  A  complete  solution  of  (2*200) 
has  thus  been  found. 

In  order  to  obtain  the  asymptotic  expansion  of 
large  |z|  we  proceed  as  before,  for  R^(z),  making  use  Of  Theorems 
I  and  VI.  From  its  definition  in  (1*57)  R2(z)  can  be  written: 

R«(z)=(lnz)z'^*- ^  v(n+v*l)z*^+z"''^  Jb  z*^+z*'^  2  b  z”  . 

n=o  xmo  ^  ns2v+l  ^ 

In  the  last  sunmatlon  put  n=2v+l-i»m»  then  from  (2*210):  ^2v+l+m^ 

=  v^(m+v+l),  so  that  for  |z|<min(la|  ,  |b|  )  we  have:  * 

f^2 ( ^)  *'^  ^  ^n^" ^ Ins)  z'^'*'”  2  V (nev+l )  z*^+z^*^  2  v.  (m+v^l )  z®- 

h=o  n=o  mio  ^ 

=z*^  Z.  b  z®+  (ins)  Z  ▼(n+vel)  Z  Vc (n+v+1)  z% 

h=o  “  n=o  n=o  ^ 


nso 


where  the  ladex  n  replaced  m  in  the  last  simiatidn  and  use  of 
(2«‘2ll)  was  made.  Thus  We  can  Write: 


2  V 

RgC  z)=z"''  2  hj^z^+ti|(  z)+U2(  z)  i  iz|<mln(  laj  ,  lb|  )  , 


(2*214) 


n=o 


where : 


(z)=(lnz)  ^(n-fvsfl)  z^^+z^"**^  %  Vg  (n4v+1)  (2*215) 

‘‘‘  a=0  n=o  ^ 

U2(z)=z'^'*‘^  X  V(n+V'fl) z“=z'^'*‘^  X  CEv^(n+v<fl)+Pv«(n+v+l)  + 


n=o 


(n*f  v+1 )  ‘fHv^  ( n-f  v+1 )  ] 


(2*216) 


Take  U2(z)  first.  For  large  lz|  its  asymptotic  expansion  is 
given  as  (2*87)  shows  for  Rj^(z)  and  for  o=v+i.  it  was  also  shown> 
immediately  after  equation  (2*114),  that: 

00  (3) 


e^«zi^/^CU^g;^^V(i2)“]^ 

=  R^(z).  Thus,  as  in  (2*87): 
v+1  4  V(*n+v+l) 

+  Re 


and  e'^^z* 


‘'*/2tl*28‘*V(-ls)“] 

n=l^ 


u, 


n^l 


.n 


rJ»z 


y*l 


^  V(*n+v+l) 


n 


*  <le 


n=l  z' 

p©  v(*n+v*l) 


R2^(z)  ,  0<4<tt 
R3(z)  ,  *tT<^<0 


rJ^z^  -  Z  — — -  +e"^®/^[GR, ( z)+HR.( z)  ]  ,  ^=0 
h=l  z“  - 

Next  take  u^(z).  It  can  he  written  as  follows: 
n“0 


(2*217) 


f  a— Q.\ 


9 


where  use  of  the  deflnltloh  (2-213)  was  also  fflade«  fhust 

U|(z)  P  (Inz)  ^  v(h‘fvtl)  2  z*' 

n=o  h=o  ^ 

in  exact  agreement  with  (2-215)*  The  part  u^(z)  in  Rg(z),  in  the 
form  shown  by  (2-218),  could  be  written  down  from  the  beginning, 
if,  in  order  to  obtain  the  second  solution  RgCz)  when  2v+l  is  an 
Integer,  we  had  followed  the  general  method  of  Probenius,  as 
applied  in  this  special  case  and  explained  In  reference  9  (pp* 
396-404) . 

The  asymptotic  expansion  for  large  |z|  of  Uj(z,o),  as  given 
in  (2-219),  can  also  be  written  down  fpllowing  (2*87),  with  a 
considered  now  as  a  variable  parameter.  Referring  also  to  (2*209) 
we  have: 

U3(z,a)  -z  -  De"  R^(z)  ,  0<(|)<tT 


^  op  v(-n+cr) 

-*'’2 — s —  +  Ce*"'/*  R,(s)  ,  -w<Ko 

h=i 

_  PP  V  ( -n+o)  /* 

'*[CR3(z)+ER^(z)]  ,  «|>=0  .  (2-220) 
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Here  a  is  a  varlatele  parameter  varying  areimd  the  root  v+1  of  the 
iadlcial  eqt^ation  about  z°d  of  equatiO!i  (I  1-50)  •  u^{z,a)  t  as 
given  in  (2-219)  for  |:z')<min(  |a|: ,  :|h|  ) ,  is  a  uniformly  oonvergent 
series  of  analytic  functions  of  &  (v(y)  was  proved  to  he  analytic 
for  all  y  and  is  an  analytic  fimction  of  o)  and  can*  therefore  * 
he  differentiated  any  numher  of  times  with  respect  to  o  (9  p* 

400).  fhenj  for  any  o,  we  let  and  obtained  the  asymptotic 

expansions  (2-220)  for  u^(z,o)  .  ^sorve  that  dependence  of _ 

vu(ZjO)  on  d  appears  now  (md  for  z  in  any  sector)  only  in  the 

^  er-  ..ft  -  ft 

series  -z  ^  v(-n4o')/z  .  v(y)  is  analytic  for  all  y  and  z°  is 
n=l 

analytic  for  all  o.  For  o^v+1  we  have:  v(-n4v4i)  =  0  (a=l,2i 

Then,  with  &  varying  in  the  vicinity  of  v+l  ^d  for 

ft  —  ft 

Sufficiently  large  |z( ,  -z  ^  v(*n+of)/z  is  a  lonlformly  convergent 

n^l 

series  of  analytic  functions  of  o*  Therefore »  differentiation 
with  respect  to  o  i^y  nisher  of  times  is  permissible,  as  before. 
Thus  we  obtain: 


lu,(z,o) 


oo  v(-n4Gr) 


Izl-po  bo*-  ^n 

Pi'bv(-n+o) 


3*.-(lnz)z^E 


nsl  z- 


•»  Z.-’sr—^  •  -"<♦<« 

n=x  z 

Finally  for  o=v+i: 


Ui(z)^%U3(z,o)3^^^^  ,~-(lnz)z 


V4i  S  v(-n+v+l) 


n^i  z 
.n 


n 


(2-221) 


^  ’bv(-n+Gr) 


_  'hv  ( -n+o) 

But:  V(-n+v+l)=a_^=  0  (n=l,2,3, . . . , )  and 


=  Vg(-n4V4i).  Therefore: 


Ut  (z)  X  ve  (“!!*►  V#i)  ^ 

^  IzHob  nsl  ^  Z 


2^60 

(2*222) 


IzV 

Sybstituting  (2»2l7)  and 


in  (2“ 


we  obtain: 


R«(z)  ^  z“^^\  z’^-z^^^E  yg(-n+v+i)ig 
2  Izkeb  n=o  ^  ^ 


1  V  (»n+v+l) 


nsl  ^ 


a=l 


=n 


R^(z)  *  6<^<ff 


XJ  .7^  ^  tr  /  -  ._V+1 

Z  ^Dv^Z  -^Z  >  Vc 

nib  nm  5 


^  V(-n+v4a:^ — 

(-n+V+l)«^  ■*Z^'*‘  ^  rrr-rr-r  jurn-BP.. 


.n 


nsl 


bH 


+ae““'^'^  R3(z)  ,  -.t:<K0 

2v 


n^o 


'rv5(-n*v+i)t  **''*^2,— T—  ♦ 

n=l  z  n=i  z 

+e  ■■■''  ■  LGR3(z)+HR4(z)1  »  ‘ 

Accor'dlng  to  (2-211)  and  (2»2lO) :  v^(-n+v4l)t^V(-n^vi*-l)sV|(»n*v^l)  = 
“^-n+2v‘fli  first  three  eunmatlons  in  the  above  e^anBlons 

can  be  combined  as  follows: 

s*'^  z t,  t  s  j”''  Jfb  zKz''*'^  ^  Saas|sai. , 

nV®  A  s'*  nfo"  nTl  s>* 

« 

Since  b  2^*’*****^  last  summation  change  the  Index 

as  follows:  n  m  2v^l-m*  The  expression  becomes: 

-v  1''.  n  «1  1E«  fm  _-v  l\  _n  ..V  S-  *’n 

z  z  -Z  -  X  2v?i-m  ®  ^^n~  ^  ^  ^ 

n=o  “  m=2v  z~  *  n=o  “  n=2v  z  - 

2Y  2v 

=  z*'^  5  b  z^»  z"^  5  b  z^  =  G  ,  l.e.  It  vanishes.  Therefore: 
h  n 

n”o  n^o 


«»  v(-n#v+i) 


R2(z) 

1  Z|-i»M 

■Sf! 

CM 

«<! 

R4(z) 

,  0  <^<TT 

^23 

R3{z) 

,  -ir<^<.0 

fSJ 

% 

R3(®) 

+  ^24  ^4^*^ 

^0 


(2-223) 


2-61 


where : 


"23 


=  Ore 


(2-224) 


The  result  is  in  cemplete  agreement  with  theory  (9  pp^  168- 
l74i  lO  pp.  72-73)  and  analogous  to  the  result  for  R|(z)i 

Determination  of  the  Coefficients  of  the  Linear  Relatians  t 
In  order  to  evaluate  Gh  and  H  we  follow  a  method  similar  to  the 
.one  Which  led  to  expresaion  (2-194)  for  D»  From  (2-111)  and _ 


(2-213)  we  have: 

Vi(y)-V|(y)=Vi(y)-v'  (y)=  ^3^(y)+FVg(y)+av^(y)«t.Hv^(y) 

As  In  (2-188)  we  can  then  obtain: 

v^(y)  V2(y)  v^(y)  v^(y)-v'(y) 


(2-225) 


M  Si 


P(y) 


v^^y+i;  Vg^y+i;  v^^y+i;  v^iy+i;-v' 
V|(y+2)  V2(y+2)  v^(y%2)  v^(y+2)-v '(y<f2) 
V|(y+3)  V2(y»^3)  v^(y+3)  Vj^(y*3)-v'' 


(2-226) 


The  procedure  that  led  from  (2-188)  to  (2-194)  shows  that  H  can 
be  expressed  as  follows: 


H  =  [vj(y+2)-v'(y+2)]p4(y-2)N^(y-2)+|^ivj^(y+3)-v'(y+3)]p^(y-l)  + 
+[v^ (y+2)  -V  ^  (y+2)  Ip^ (y-l)|  N^(y-l) -^v^  (y)  -v '  (y)+2a[vj^ (y+l)  - 
-v'(y+l)  ]jN^(y)-[vj(y-^l)-vMy+l)  ]N^(y+l)  ,  (2-227) 

Differentiate  (2-194)  with  respect  to  y: 


0  =  f' (y+2)p^(y-2)N^(y»2)+V^(y+3)p^(y»l)N^(y-l)«^v'(y+2)p5(y»l)  ^ 
•N^(y-l)-|:v'  (y)+2av'(y+l)  ]N^(y)«v'  (y4l)N^(y+l)+v(y+2)|*Cp^(y-2)» 
^N^(y-2)  ]+v(y+3)|^[p^(y-l)N^(y-l)  M(y^2)|^[p^(y-l)N^(y-in- 

-Cv(y)+2ay(y*l)]CdN^(y)/dy3-v(y+l)[dK^(y+l)/dy3  ,  (2-228) 


i«62 

Add  and  (2-228)  : 

H  =  Vj^(y+2)p2^(y*2)N^(y-2)'^Cvj^{yH‘3)p4(y-l)+v^(y+2)pj(y-l)i:N^(y-l)- 
»Cvj^(y)+2avj^(y^‘i)  3N^(y)-Vj(y‘^i)N^(y+l)+v(y+2)|yCp^(y»2)  • 

‘N^(y“8)  ]‘*‘V(y+3)|“Ep2^(y-i)N^(y-l)  ]+v(y+2)|^ipj(y-l)N^(y-i)  ]- 
Cv(y)#2ayCy+l)  ][dN2^(y)/djJ-v(y+i)CdN^(,y4i)/dy]  *  ( 2-229) 

Kitting  y  =  n+vH-l,  whene  n  any  integer  *  and  Gbserving  that: 

v^(n+v+l)  =  B^=  ^n+2v+l  ^  ^(n+v+l)  =  a^^  , 

we  obtain: 

H  =  in+3P4(n+v»i)N^(n+v-l)+Clj^^3p4(n+v)+B^^2p^(n‘fv)  ]N^(n+v)-[l^+ 

*2ai^^j}N^(n+v+i)-Bn^lN^{n+v+2)+an^2[p|^(n-i‘v»i)N^(n+v-l)  + 

+p^(n+v^l)Nj^(n^v«i)  ]+a^^^Cp^(n+v)N^(n+v)+p^(n+v)N^(n*fv)  ]+ 
♦aj^^2tp3(n+v)N^(tt4v)+p3(n+v)N^(n+v)i-[aj^+2a  aj^^^iN|(n^v+i)- 
♦^an^.^N|(n+v+2)  , 

or: 

.p3(n+v)+a^_i_3P^(n+v)+a^^2P3(n+v)  ]N^(n+v)  *[B^+2aB^^^]N^(n+v+l) - 
"®n+1^4  ^  n+v+2 )  +a^^2P4  ( n+v*l )  N  ( n+v*l )  +[  a^^^p^  ( n^-v )  ( n+v )  |  • 

•Nj[(n+v)«[a^*2a  aj^^-L]N|(n+v+l)^ajj^3^Nj^(n+v+2)  ,  (2*230) 

Exactly  the  same  expression  gives  0»  if  N3(y)  is  substituted  in 
place  of  N^(y) . 

An  expression  for  N'(y),  at  least  asymptotic  for  large  |y) 

In  the  right  half  y*plane,  can  be  obtalhed  by  direct  differentia- 


tlon  of  (2«16a) -(2-169)  o?  (2-178) -(2-179)  (17  pp*  434-461 
PP»  457*459).  We  obtain*  eorreepondlngly ! 

<(y)  =  ttiw/2.*(y.v.4).t(y)}»^(y).  ilHl 


4' 

3 


2  r(y+v-»»3?lG/2) 


.i»(y.v.3?u/a).  , 

y+v+3lie/2 


dgt(y»v-f5lic/2) 


(y+v+3+lc/2)  (y'fv+4?lc/2) 


. . ] 


(2-231) 


N4(y)  =  C+iff/2+'jf(y+v+4)+iKy+3-v)]N.  (y)  + 
3  3 

r,/  /.  dT^Ky+7+ic/2) 

y+6?ic/2 


(+1)^  r(y+v+4)  r(y+3»v) 

r(y+6?lc/2) 


4g^(y*8lic/2) 
(y+6lic/2) ^y+7+ic/2) 


* . 3 


(2*232) 


where  t(z)  =f'(z)/r(2)  (17  pp*  241-267)4  n%)  satisfies  the 
relations: 


t(z+l)  =  t(z)+  4 
z 


'KD  =  -Y  -  tier's  constant 


(2-233) 


It  be  checked  easily  that  for  real  a,  b:  G  =  H.  We  have 
expressed  H  (G)  in  terms  of  N^(y) ,  Nj^(y)  (N^,  Np  only,  just  as 

D  (C)  was  expressed  in  terms  of  N^(y)  (N^). 

Concerning  the  computation  of  C,  D,  G,  H  through  equations 
(2-195)  and  ( 2-230)  we  observe  that  n  can  be  given  values  large 
enough  so  that  the  factorial  series  for  N(y)  and  N'(y)  are  easily 
and  quickly  computed.  For  such  n,  r(x)  and  ^(x)  can  be  evaluated 
with  the  use  of  their  well-known  asymptotic  expansions  t 


2*64 


Even  for  x=2  tlae  seGOftd  series  yields  ^(2)  with  an  accuraGy  of 
6  deGimals,  while  the  same  is  true  for  r(x)  and  x=3» 

Notioe,  however j  that  n  oan  not  be  given  very  large  values 
in  (2-195) »  (2~230)  *  Since  ^d  B^^'s  (the  coeffiGients  of  the 
power  series  in  x  of  R^(x)  and  R^Cx)  )  are  involved  in  these  rela» 
tions  and  the  aoeuracy  of  their  evaluation  diminishes  as  n  inorea* 
aes,  there  is  a  limitation  to  the  values  of  n  that  can  be  used* 

It  was  also  observed  that  for  large  n,  the  summation  of  the  terms 
in  the  right  half  sides  of  (a^^l^S)  and  (2»230)  destroyed  the 
accuracy  rapidly  by  eliminating  the  first  significant  decimals  of 
the  individual  terms.  In  each  particular  case  there  is  an  optimum 
range  of  values  of  n  for  which  (2-195)  and  (2-230),  with  a  given 
accuracy  of  computation,  yield  the  most  accurate  results.  To  make 
sure  that  the  values  of  C,  D,  0,  H  are  the  correct  ones,  one 
should  use  (2=195)  and  (2=230)  for  about  10  values  of  n  and 
compare  how  well  the  10  values  of  these  coefficients  agree. 

As  an  Indication,  we  give  a  few  results  obtained,  with 
8=declmal  accuracy  machine  computations,  In  Case  I,  Chapter  3, 

PART  I.  In  this  case  a^l2,  b=lO,  C-D,  0=1,  Rj(x)=R^(x).  iesides 
p  and  H  ,  the  values  of  BgC??)  for  x=i2  and  x=l4  are  given  for  the 
first  five  functions:  v=i,3,3,7,9.  R^(x)  is  not  given,  since  it 
is  Included  In  (1-80),  defining  R2(x),  x=i2  and  x=14  fall  In  the 
overlapping  region  between  the  convergent  series  (1-80)  for  RgC^) 


....  ,  -'i.  .  , 
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and  tne  asjraptotiG  sanies  (2-223)  J  R2(x)=e“^®'^^[GRj(x)+iBR^(x)  ] 

yielding  values  of  RgCx)  thnoi^ii  §#  H  and  tne  asymptotle  series 
(1-24)  for  R3(x)  and  R^(x).  Both  valuefl  of  Rg(x)  are  given  for 
comparison. 

V  D  M 

1  4.14617-14.06427  -252*508-1265.563 

3  -I33.l93««‘ill7*365  _  a98*741+il306*218  _ 

"5  (1.24526-10*965800) *10^  (-0.875955-14*30989)^10^ 

7  (-2.32660+11.5730) 40^  (-2.05026+13*83348)^10^ 

9  (7.2577-14.2497) *10®  (10.0441-17*42212)^10® 


V 

X 

B.^(x)  from  (1-80) 

1 

12 

14 

38.0483 

114.081 

12 

-327.000 

3 

14 

-401.884 

12 

1.46444.10^ 

5 

14 

7.44758*10® 

12 

-1.95279.10® 

7 

14 

-2*95473*10® 

12 

4.80712*10® 

9 

14 

-3*75249*10® 

from  (2-2223) 


38*1064 

114.082 


-326*902 

-401.875 

1*46271 *10^ 
7.4479540^ 

-1. 95571*10® 
-2*95481*10® 

4.80351*10® 

-3*75230*10® 


For  x=14*  falling  roughly  In  the  middle  of  the  overlapping  region 
In  this  case,  the  agreraent  is  good  up  to  five  significant 
decimals. 


Fig  2-2  Integration  paths  (3  ,(4  for  V3(y)  ,V4(y) 


Fig.  2-5  integration  paths  13,(4  Fig.  2-6  Integration  paths  (3, 

for  N3(y)  ,N4(yX  for  N3(y),N4(y)  . 
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CHAPTlR  3 

imAmS  AND  aSNSRALIZAflONS 


ffee  analysis  given  in  the  previous  chapters  led  to  the 

complete  solution  of  the  proiblem  for  the  special  forms 

of  the  stratification  function^  The  complexity  of  the  prohlem 
depends  exclusively  on  the  form  of  cip(x)i  since  it  is  this  form 
which  determines  the  numher  and  nature  of  the  singularities  of  the 
radial  equation  (l  1*42).  There  are  three  other  forms  of  ^(x) 

for  which  the  method  used  for  (|>(x)=  can  be  readily  applied. 

They  are  given  below,  together  with  the  differential  equation 
into  which  (I  1*42)  reduces,  respectively: 


s  0  (J-l) 

(x+a)(x+b)  -  °  x'^ 

(^(x)=  ;  puts  A  z  ,  (j>(z)  =  m 

♦  —Jd— In' (*)♦[— as—  .  IiI^]R(z)=  0  (3-a) 

--  (z+a)(z+b)  4z(z<fb)  4z- 


:  p«t:  A  ,  , 

x~+b 

R  (s)*[fe  +  - - - - -jR 

(z+a)(z+b) 


<tH),  [|jt] 

(z+a)^  v(v*l) 
4z  ( z+b )  ^  4z- 


z)=  0  (3-3) 


All  the  above  forms  of  represent  stratifications  similar  to 

(x):=  Shown  In  figure  (1-1) ,  PART  I.  The  essential  point 


5^a 

is  that  the  radial  differential  equation  for  TM  waves  ohtained  in 
all  cases i  has  the  same  numher  and  nature  of  singularities  as 
(I  l«SO) .  More  speGlfically i  three  regular  singularities  at  the 
finite  X  (or  z)  plane  (one  at  x=0»  or  z=0)  and  an  irregular 
singularity  of  the  first  rank  at  6o  *  in  the  last  two  cases  the 
solutions  around  z-oo  hecome  suhnormal  in  nature  (9  PP*  4l7“4t8 
16S“171,  10  pp«  63*64),  hut  the  modifications  required  are 

elementary.  She  correspohding  difference  equation  (il  2*5 )  can _ 

again  he  solved  along  identical  lines.  Ford's  method  can  he 
applied  equally  well  to  these  eases  to  provide  the  precise 
asymptotic  expansions  required  (  7  Ghapt.  vill) *  in  addition  to 
Ford's  fheorems  l  and  VI,  used  in  the  preceding  analysis,  we  may 
have,  in  the  last  two  cases,  to  make  use  of  a  number  of  other 
similar  theorems  contained  in  his  hook,  in  the  last  chapter, 
reference  7,  relative  examples  are  included,  showing  that  the 
method  applies  without  essential  modifications.  It  must  he  pointed 
out  that  the  analysis  presented  in  the  previous  chapters,  has 
modified  and  generalized  Ford's  method  in  two  directions*  It  has 
provided  general  expressions  for  the  coefficients  G,  P  and  G,  H 
depending  only  on  one  solution  of  the  adjoint  difference  equation 
(II  2*137)*  And,  mainly,  it  has  dealt  successfully  with  the  case 
of  integral  values  for  the  difference  of  exponents  (i*e.  when 
2v+l  is  equal  to  a  positive  inteder) ,  almost  always  present, 
directly  or  indirectly,  in  all  physical  prohlems.  Both  these 
generalizations  are  applicahle  in  the  last  two  cases  (3*2)  and 
(3*3)* 

The  modification  of  Ford's  method,  to  which  we  Just  referred, 
consists  of  introducing,  in  place  of  x,  the  variable  y=x+cr  and 
treat  the  difference  equations  (II  2»5)  and  (li  2*137)  in  terns 
of  y.  The  parameter  c  is  always  present  in  Ford's  work  (called  h) 
and  makes  the  extension  of  the  analysis  to  integral  values  for 
2v+l  almost  impossihle.  By  introducing  y  we  eliminate  this  para* 
meter,  without  rendering  the  method  Inapplicahle;  at  the  same  time 
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we  are  able  to  extend  It  easily  to  the  case  when  2v+i  takes  on 
integral  values ^ 

ffae  question  that  finally  arises  is  whether  the  analysis 


developed  in  the  preGeding  chapters  can  he  applied  iuccessfully 
to  more  general  types  of  stratification;  for  in8tance»  to 
varying  in  a  manner  similar  to  figure  (3^*1)  • 


Figure  3*i»  More  general  forms  for  c|>(x) 


Almost  always,  (except  in  trivial  cases),  the  answer  depends 
on  the  nature  and  number  of  singi^aritles  that  are  introduced  in 
equation  (I  l«-42).  The  point  is  .illustrated  hy  mentioning  certain 
cases  where  the  method  fails.  If  Irregular  singularities  are 
introduced  in  the  finite  x  (or  s)  plane,  then,  in  spite  of  the 
fact  that  power  series  solutions  can  still  he  found.  Ford's 
theory  for  obtaining  their  asymptotic  expansions  is  no  longer 
applicable  (  7  Cbapt.  Vill).  The  difference  equation  (11  2^5)  is 
no  longer  normal  and  its  solutions  are  of  such  complexity,  that 
no  corresponding  theorems  (like  I,  VI,  or  the  rest  of  Ford's 
theorems)  exist,  which  can  he  applied  to  these  solutions  to  yield 
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tke  asymptotic  expansions  Or  R^(x) *  ^^(x).  A  case  like  this 
arises,  if  <^(x)  =  reason  why  it 


was  not  included  in  the  previous  list  (3-1)  to  (3-3)*  It  is 
easily  seen  that  x=»a  ^d  x=*fe  are  irregular  singular  points  of 
(1  l-4g). 


Another  case,  where  the  method  fails,  arises,  if  the  singu¬ 
larity  at  x=6o  (or  at  least  z^o©,  if  we  must  resort  to  the  _ 

change  of  variable  x=x^^j  ),  has  a  rank  higher  than  1,  Normal 
asymptotic  series  for  i^(x),  R2|^(x)  can  still  be  obtained,  in 
general,  but  Ford's  theory  again  fails  (  7  pp*  339-341),  owing  to 
the  complexity  of  the  solutions  of  the  difference  equation  or  to 
the  impossibility  of  even  solving  it^ 


Another  ease  arises,  if  the  number  of  finite  regular  singu* 
larltles  is  large  or  infinite.  For  example,  Cf>(x)  s  l+ae"^^, 

(b  >  6) ;  the  coefficient  (^'(x)/<^(x)  of  R'(x)  in  (I  1-42)  becomest 


-abe 

l^ae 


-bx 

-bx 


,ab^ 


a*6 


bx 


and  introduces  an  infinite  number  of  regular 


singularities. 


at  the  zero©  of  e^^^a, 


in  the  complex  x-plane. 


Furthermore,  in  this  case,  the  irregular  point  xsqo  can  not  be 


assigned  a  finite  rank* 


In  practice  such  cases  can  be  treated  either  numerically  or 
by  approximating  the  function  Ci>(x)  by  more  simple  fvmctions*  The 
numerical  results  in  Chapter  3*  3P^T  1,  showed  that  this  approxi*- 
mation  is  valid  and  permissible.  For  sharper  variations  of  a  more 
complicated  nature,  one  may  divide  the  interval  O^x^op  into  a 
finite  number  of  shorter  intervals*  The  problem  then  is  somewhat 


similar  to  stratification  by  layers,  requiring  additional  matching 
processes  at  each  spherical  boundary  separating  regions  of 
different  functional  representation  for  ^(x). 

If  the  stratification  terminates  at  a  finite  distance  x. 


series  solutions  may  be  suffieient,  as  it  has  already  been  pointed 
out;  anyway,  the  problem  can  be  classified  as  a  special  case  of 
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Bts^atification  by  layers. 

It  must  also  be  pointed  Out  tbat  there  are  problems  where 
only  fi  waves  are  involved^  The  radial  equation  (I  1-43)  for 
these  waves  Is  more  simple  than  equation  (1  1-42)  for  TM  waves# 

In  certain  cases ^  the  former  equation  can  be  solved^  while  the 
latter  can  not#  An  example  is  mentioned  In  reference  but  it 
refers  to  a  finite  interval. 

With  these  remarics  In  mind  we  can  answer  the  question  raised 
Cfor  problems  requiring  solution  In  the  whole  Interval  0  4x4^)f 
as  follows:  The  preceding  analysis  can  be  applied  to  more  general 
stratification  functions  ^(x)  as  long  as: 

1)  No  finite  Irregular  singularities  are  introduced  in  equation 
(I  1-42)* 

2)  The  rank  of  the  irregular  singularity  at  x  or  z=«o  does  not 
become  higher  than  1;  in  general,  this  requirement  is  more 
easily  satisfied  than  the  preceding  one,  since  In  all 
physical  problems,  (|(x)  must  reduce  to  1,  or  a  finite 
constant,  at  x  or  z=o6* 

3)  The  number  of  regular  singular  points  introduced  in  the 
finite  X  (or  z)  plane  Is  reasonably  small  and  permits  the 

use  of  a  bilinear  change  of  variable  x  (or  z)  =  y^f  »  which 

maps  the  interval  of  interest  in  the  x#plane  (in  real  cases 
the  real  x«axls)  Into  a  finite  circle  arowid  t^-p/a  In  the 
t-plaae,  placing  all  the  other  singularities  of  the  equation 
outside  this  circle.  This  last  restriction  is  not  as  funda¬ 
mental  as  the  first  two,  but,  from  the  computational  point 
of  view,  it  marks  the  difference  between  the  possibility  of 
an  analytical  or  completely  numerical  solution  to  the 
problem. 
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CbwM  K<M|  ValvaraUy 

Eltitrltal  Ea«laa*rlk|  DapaMataU 

Ttteaa,  Talvaa 

EapablU  al  Cklaa 

Alla.  Pralarrar  Ckap'NU  CA*« 


Mr.  D.  $■  laaat 
OapartniaM  al  Mf.haaiaitr* 
Uaiv.  CaltoaaalNa.  AullarAaM 
Raala.  Mailariablr*,  Ekilaal 

Pralaaaar  Paal  Barul  WUa 
Oaaka  City  Uaivariliy 
Oapi.  al  Caytatarlap  AcUaaar 
II  Niaki  Oplmaibl  KUaba 
Oraka,  lapaa 

OaMl4  C  Mlaaaa 
Oapi.  al  CU<iri«al  Baar. 
UalvPifli*  al  AritPNk 
MrM  ti,  ArliaM 


Maala*!  Aaaaarib  lauliala 

Wrala  Park,  CaUlaraia 
AMa:  Maty  Laa  PiaUi 

Or.  C  H.  P»pa* 

Dipt  al  Clatirlial  EaaiaairlaA 
Cablaraia  laaiiialk  pl  TarkkkiPff 
PaiaAaM,  Caltlaraia 


•taaNt!,  CaUlaraia 
Alta:  DacMkiM  Library 
Applil!  CIltirPM** 


Protkikkr  IfrpM  E.  UaBar 

01*.  kl  Cla«<rt*al  SaAlkiarlaf 
Uaivariliy  *1  CallMrata 
Earbalay  4,  CalllarM 

Pralairar  CharU*  lUHal 
DaparUkaat  *i  Pkyaic* 

UMvarrUy  *1  CaillarMa 
Earbalay  4,  Calltorato 

•arui*  library 
EraaAali  Ualvarally 
Vabkaa,  MapaarkaaaNp 

PraNaaar  H.  O.  Eaakar 
Ackaai  al  ElaairUal  EafiaaarbM 
Caraall  UMvarrlly 
bkata,  Nav  Tark 

Library,  CaUapa  al  O^kaartkA 
Uki>«r*uv  HaitUf  Library 
Ukivariliy  NalfMa 
NkM  Tark  Ualvarrilp 
N*m  Tark  LI.  Nk*  Tark 

C  A.  Chipiaia.  Likrarbta 
Eaaaaalaar  PainaahM  laaNNNk 
Aktai  Eala*  Hall 
Trpy,  Na*  Varh 


Oaparukiai  al  ElatlrUal  E^aakriiM  Eaaart  Plaaiay 

Till  Uai***iil|t  Oaparlatau  al  CaiiavartaA 

btvHarra.  CaMttli**  C***  lariUali  kl  faakaili^ 


Aaiiiiiaa  Labarilary 
Lakai  Hapkiar  Ukivarrily 
ILI I  N.  Pkal  Mraal 

AilWMri  1.  Ukryla^ 

Ptirttar,  Liarkla  L*k*iklary 
Mill  lariilvla  a*  Titkaalapy 
A*4la*4.  UiiiiibaMMi 

Mr  Lak*  H*«iM 
DaramaM  Aavtn 
Aaarartk  Lab  *1  Eiartrarutr 
M»*»  lariiivl*  *1  Titkaalapy 
Caink*i4pi  >4,  Mariatkaialli 

P*vl*>***  *  Va*  Mippal 
Man  IkriiiM*  *1  Taikaalaay 
lab-  la*  taialallaa  Aiiaartk 
Caa<brl4p*  >4,  Mkiiarbaraui 

library,  kawa  A  114 
Liarala  Lakarktaiy 
P  O  A*>  ti 

laaiap'aa  IL.  t4*>i*«lt*i*ii* 


Dap*.  alClaalrltii  Ei^aaarlap 
VA**  Uiltlali  *1  Tattaaiapy 
UMvariiiy  Ctrila 

CU**taa4  4,  Okt* 

All*.  1  Aaaly.  HaaA 

Or.  e  L  Pilkaakark 
Aalall*  MaatarUI  IkailiM* 

Calukbak.  OM* 

Alta:  ElatirKal  EaylMaHkp  DlvtiUk 
Likrariaa 

Baptaa**iap  Library 
Ar*«k  Uaivariliy 
PravtAaa**,  AbaA*  lala^ 

Praia****  A  *  liratiaa 
Dap*.  *1  CirtitirkI  Capiaaana# 

Uaivariliy  *1  Tiiai 

Avail*  11.  Tiiai 

Ur  VilUaai  Way 
Eaaaaitk  Librkrlk* 

T*»k*r  Uiirvaakai*  Ckiy 
1>44  Cakvaap*  Aaalavar! 

Hallrtrav*  II.  Calilaral* 


Tbaary  aii  Aaalyvl*  Oapartmaal 
Will—  Aua  Lakaraiariai 
Vaiwrviiy  *1  Michryaa 
VtllBv  Ava  Airpari 
Tpailaati.  Uithlpaa 

•OUO  •TATE  ONLY 


Pivantapaaiic  Piiilaa 

•all!  fiK*  Piviiiaa 
Navil  Aiiairck  Lakaraiary 
Vaibiapia*  II.  O  C 
Ana.  C*4a  4411 

Or.  Aaiarar  •—■***•,  Lr 
Otlaaat*  Mtiirlali 
B***irc4  Lakarilary 
WatiriawB  AiviMl 
Waiaria— ,  UaiiathwaMti 


Clatiraaixt  Aaiikrrk  lAbariiary 
Oiviiiaa  al  Elaclrital  Baplaaatlat 
Uaivariliy  *1  Calliarala 
Airkvlir  4,  Cilltaral* 


Lalwir  Hlyklai  Uaivariliy 
>41*  in!  Charlii  Stint 
Vkliik**!  Hill 
Aaltiinari  II,  Mirylaa! 
All*  Ur  L  O  Alima* 


PkyiKi  DipirlmiM 

Amharll  Callipi 
Amkiril,  Miiiathaiitii 


Pralniar  A  C  Harharp 
Daparimrai  al  PkyiKi 
Wiiblaptaa  Uaivariliy 

M.  LavIi,  Uiiiawri 

uirtavravi  Aiiaartb  IkatlNA* 
PalyiachM*  laHiivl*  *4  Eraablyk 
II  Lakaiaa  Strait 
Araaklyn,  Ni«  Tark 

Alta:  Librartaa 

Dr  SIAiiay  Sbaptra 

Arlbvr  0  Lilli*.  U*. 

II  Atari  Park 

CkJtikrUp*  40,  Uliaaibaian* 


Mr  WiUiiaaH  Pram 
Eva*  KMpM  Carparatlan 
104  A  Mraal 

NiaSlwra,  M****ck«**M* 


Or  W  U  Wklib 

Ball  Talipbaa*  Ub*  ,  lac 

Uvrriy  Hill,  Mav  Lanay 


Dvparlmail  al  Phyii 
UaivitiHv  *1  Ukanai 
UlaMipaiii,  Mlaaai 


ISU  Will**  lAbaiaiail** 
III  Will  lllik  Mrial 
Ha*  Tark  IT,  Ni*  Tark 


ELCCTAOMAONETIC  AADIATION  ONLT 


Palyiixkax  ininiuli  *1  Draaklya 

ll  Lvkn.ai  liiwl 

AiaaUrn,  Ni«  Tark 


Pralaiiai  O  E  N  k| 
b  O  •*■  144 
Ailmir,  Nav  Lariiy 


■  patllli 


